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ABSTRACT 


In  the  past,  human  body  models  have  been  developed  by  assuming  simple 
geometric  shapes  for  the  components  of  the  human,  body  in  order  to  predict  the 
mass  properties,  i.e.  mass,  center  of  mass,  and  inertia  tensor.  In  this  study 
a  new  personalized  method  of  predicting  component  mass  properties  is  developed, 
based  on  experimental  data  for  whole-body  mass  properties  for  three  body 
positions  where  only  one  component  is  moved.  Also  needed,  however,  are  esti¬ 
mates  for  the  mass  and  sum  of  the  moments  of  inertia  of  each  component. 

Additionally,  the  contribution  of  outer  garments  or  a  suit  to  the  mass 
properties  can  be  included  in  the  prediction  if  it  is  assumed  that  the  compo¬ 
nents  of  the  outer  covering  are  rigid  bodies  that  move  negligibly  relative  to 
corresponding  body  components. 


iv 


TABLE  OF  CONTENTS 


Page 


Acknowledgments  .  iii 

Abstract . iv 

List  of  Figures . vii 

Chapter  1  Introduction . .  .  .  . .  1 

Purpose . 1 

Background  .  2 

Study  Structure . .  .  .  3 

Chapter  2  Determination  of  Whole-Body  Mass  Properties  from  Component 

Mass  Properties  with  Respect  to  Component-Fixed  Axes  ...  5 

Notation  .  5 

Analysis . 9 

Chapter  3  The  Hanavan  Model  . .  13 

Chapter  4  Determination  of  Component  Mass  Properties  about  Component- 

Fixed  Axes . 16 

Overview . 16 

Component  Centers  of  Mass.  . . 17 

Component  Inertia  Tensors.  .  .  .  25 

Experimental  Determination  of  Whole-Body  Mass 

Properties . 36 

Consideration  of  Outer  Garments  or  Suit . .  37 

Chapter  5  Explanation  of  the  Infinity  of  Solutions  for  Component 

Inertia  Tensors . 39 

Chapter  6  Suggested  Methods  for  Estimation  of  Component  Masses  ...  43 

Regression  Equations  .  43 

Water  Immersion . .  45 

Inclusion  of  Outer  Garments  or  Suit  in  Component 

Masses . .  46 


v 


Page 


Chapter  7  Suggested  Methods  for  Estimation  of  the  Sums  of 

Component  Moments  of  Inertia  .  .  47 

Hanavan  Model  . .  47 

Casts  of  Subject's  Components .  48 

Inclusion  of  Outer  Garments  or  Suit  in  s .  48 

J 

Conclusion . .  .  .  . .  49 

Chapter  8  Limitations . 50 

Invalidity  of  Assumptions  . 50 

Possible  Errors  in  Whole-Body  Experimental  Data  ...  52 

Accumulated  Errors  in  Central  Component  .  52 

Joint  Constraints .  53 

Suggested  Sensitivity  Analyses  .  53 

Chapter  9  Conclusion  .  55 

Appendix  A  Computer  PROGRAM  WHOLE . • .  58 

Appendix  B  Computer  PROGRAM  PART  .  62 

Bib!  iography . 71 


Vita 


LIST  OF  FIGURES 


Figure  Page 

1  Vectors  and  Coordinate  Systems  .  8 

2  Rotation  Angles . ■ . . .  .  11 

3  The  Hanavan  Model . 14 

4  Flow  Chart  of  the  Determination  of  Component  Mass  Properties  .  .  57 


vi  i 


CHAPTER  1 


Introduction 


Purpose 

The  motion  of  a  dynamic  system  depends  on  the  mass  properties  of  the 
system,  i.e.  the  mass,  center  of  mass,  and  inertia  tensor.  If  a  human  being 
is  a  part  of  the  system,  then  the  system  properties  depend  on  his  or  her 
contribution.  The  human's  contribution  to  mass  properties  depend  on  body 
position,  and  hence  become  a  function  of  time  if  the  human  is  in  motion. 

For  many  dynamic  systems  the  human's  contribution  may  be  negligible, 
as  for  a  large  airplane,  where  the  pilot's  mass  is  small  compared  with  the 
mass  of  the  system.  However,  there  are  systems  where  this  is  not  the  case, 
for  example  an  astronaut  with  a  pressure  suit  and  backpack.  For  some  dynamic 
systems,  the  human  may  constitute  the  entire  system,  e.g.  a  swimmer,  a  diver, 
a  gymnast.  Where  the  human's  contribution  to  the  system  is  significant,  a  pre¬ 
diction  of  the  system  dynamics  depends  upon  knowledge  of  the  time  history  of 
the  person's  mass  properties.  While  a  human's  mass  is  directly  measurable 
(using  scales),  this  is  not  the  case  for  the  other  mass  properties. 

The  purpose  of  this  study  is  to  develop  a  method  for  predicting  the 
center  of  mass  and  the  inertia  tensor  for  any  human  being  in  any  given  posi¬ 
tion,  with  or  without  outer  garments  or  suits.  Thus,  in  this  study,  methods 
are  presented  for  obtaining  the  time  history  of  the  mass  properties  of  any 
particular  human  being  directly  from  the  time  history  of  the  body  position. 

The  means  presented  for  predicting  body  center  of  mass  and  inertia  tensor 
shall  depend  upon  the  following  assumptions: 
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(1 )  The  human  body  is  adequately  modeled  by  a  system  of  rigid  bodies 
connected  by  joints. 

(2)  There  exist  valid  method  for  estimating  the  mass  of  each  body  com¬ 
ponent.  A  review  of  possible  methods  is  given  in  Chapter  6. 

(3)  There  exist  valid  method  for  estimating  the  sum  of  the  three  moments 
of  inertia  for  each  component  about  component-fixed  axes.  Possible  methods 
are  discussed  in  Chapter  7. 

Assumption  (1)  above  is  equivalent  to  assuming  that  body  components  do 
not  deform  appreciably  for  different  body  positions.  This  assumption  is 
questionable  for  the  torso,  whose  curvable,  twi stable  backbone  can  make  seem¬ 
ingly  significant  changes  in  torso  shape.  For  a  suited,  backpacked  astronaut 
this  assumption  seems  more  reasonable  because  of  the  restricted  mobility  im¬ 
posed  by  currently  used  equipment. 

Assumption  (2)  arises  as  a  consequence  of  the  fact  that  there  appears 
to  be  no  nondestructive  means  for  experimentally  determining  the  mass  of  the 
individual  body  components.  There  seems  no  way  to  directly  weigh  the 
individual  body  components  of  a  living  human  because  of  the  unknown  forces  at 
the  joint(s)  connecting  each  component  to  the  rest  of  the  body. 

Assumption  (3)  is  necessary  to  provide  a  needed  additional  independent 
equation  relating  the  component  moments  of  inertia.  The  necessity  for  the 
introduction  of  this  equation  is  established  in  Chapter  5. 


Background 


Several  researchers  have  attempted  to  obtain  experimental  values  of  the 
mass  properties  for  parts  of  or  entire  humans.  Much  of  the  early  work  in  this 
area  consisted  or  obtaining  centers  of  gravity  and  moments  of  inertia  for 
parts  of  dissected  frozen  cadavers.  Thus  Braune  and  Fischer  obtained  these 


values  for  a  sample  of  three  (Refs*  2,  3),  Fischer  later  added  another  one 
(Ref.  8),  and  Dempster  added  eight  more  (Ref.  6).  Barter  used  the  data  thus 
obtained  to  formulate  a  set  of  regression  equations  estimating  the  weight  of 
each  cadaver  component  as  a  function  of  total  weight  (Ref.  1). 

Later  experimental  work  centered  on  mass  properties  for  living  humans. 
Swearingen  determined  the  centers  of  gravity  for  5  living  subjects  in  67 
body  positions  (Ref.  16).  King  investigated  the  locus  of  the  center  of  gravity 
for  various  body  positions  (Ref.  12).  Santschi ,  DuBois,  and  Omoto  determined 
the  center  of  gravity  and  moments  of  inertia  for  66  subjects  in  various  body 
positions  by  using  a  complex  pendulum  (Ref.  15). 

Fowler  mathematically  showed  how  all  the  mass  properties  of  a  spacecraft 
could  be  experimentally  determined  in  flight  by  using  three  independent, 
known  thrusts  (Ref.  9).  This  method  is  applicable  to  a  human  subject,  and 
is  the  only  available  experimental  means  of  determining  products  of  inertia. 

More  recently,  attempts  have  been  made  to  model  the  human  in  order  to 
predict  mass  properties  as  a  function  of  body  position.  Hanavan  (Ref.  11) 
assumed  a  simple  geometric  shape  for  each  body  component  and  derived  formulas 
for  inertia  moments  about  component-fixed  axes.  DuBois  and  others  improved 
Hanavan' s  model  and  added  a  space  suit  model  consisting  of  hollow,  simply 
shaped  "shells"  (Ref.  7).  Tieber  and  Li ndemuth  (Ref.  17)  made  further  modi¬ 
fication  of  Hanavan's  model  and  also  modeled  a  Gemini  space  suit.  They  used 
an  improved  set  of  regression  equations  developed  by  Clauser  and  McConville 
(Ref.  4)  to  obtain  estimates  for  component  masses . 

Study  Structure 

This  study  first  shows  how  whole-body  mass  properties  can  be  determined 
from  component  mass  properties.  In  this  procedure  component- fixed  axes  are 
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employed  (Chapter  2).  Then  two  methods  for  obtaining  these  comp  orient  pro¬ 
perties  are  presented.  The  first  of  these  is  a  model  formulated  by  Ernest 
P.  Hanavan  (Ref.  11),  which  assumes  simple  geometric  shapes  for  the  components 
(Chapter  3).  The  second  method  is  a  more  general  approach,  making  no 
assumption  about  component  shapes  except  that  they  are  rigid  (Chapter  4). 

The  necessity  for  the  introduction  of  an  additional  equation  relating  the 
component  moments  in  order  to  implement  this  latter  method  is  then  established 
(Chapter  5).  This  second  method  depends  upon  experimental  data  for  whole- 
body  mass  properties  for  selected,  fixed  positions  (Chapter  4,  page  36), 
estimates  of  component  masses  (Chapter  6),  and  estimates  of  sums  of  component 
moments  of  inertia  (Chapter  7). 

A  discussion  of  limitations  of  the  second  method  and  recommendations 
for  further  study  are  presented  (Chapter  8),  and  finally  a  conclusion  is 
given  (Chapter  9). 


CHAPTER  2 


Determination  of  Whole-Body  Mass  Properties  From 
Component  Mass  Properties  with  Respect  to 
Component-Fixed  Axes- 


Notation 

The  notation  for  this  chapter  and  for  the  remainder  of  the  study  is 
as  follows: 

Scalar  Notation 

i , j  -  Subscripts  denoting  the  quantity  subscripted  is  for 
component  number  i  or  j,  respectively . 

mD  -  the  total  body  mass. 

D 

m..  -  the  mass  of  component  i  . 

n  -  the  total  number  of  body  components. 

s.  -  the  sum  of  the  moments  of  inertia  of  component  j 

J 

about  an  axis  system  centered  at  j's  center  of  mass. 

X,  Y)  p  -  a  set  of  rotation  angles  between  the  axes  of  C.!  -  xyz 
and  c.j  -  xyz  (these  angles  are  more  explicitly  defined 
later). 

I,  II,  III,  ...,  N  -  superscripts  denoting  the  quantity  superscipted 
is  for  a  body  position  denoted  by  the  Roman  Numeral  (  N 
denotes  the  largest  body  position  number). 

)  -  an  operator  denoting  the  difference  between  the  value 
of  the  quantity  within  the  parentheses  evaluated  for 
position  I  and  the  value  for  position  II  (any  Roman  numerals 
can  be  used) . 
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Point  and  Coordinate  System  Notation 


A 

A-xyz 

CB 

Ci 

Ci-xyz 
C!  -xyz 


J-xyz 


J-  -xyz 


an  arbitrarily  located  (but  known)  point. 

an  arbitrarily  oriented  Cartesian  coordinate  system 

centered  at  A  . 

the  center  of  mass  of  the  body. 

the  center  of  mass  of  component  i  . 

a  Cartesian  coordinate  system  centered  at  C^  with 

axes  respectively  parallel  to  the  axes  of  A-xyz  . 

a  component-fixed  Cartesian  coordinate  system  centered 

at  C.j . 

the  joint  (pivot  point)  at  the  end  of  component  i 
closer  to  the  body  torso  (e.g.,  the  elbow  is  the  fore¬ 
arm's  J.j  ). 

a  Cartesian  coordinate  system  centered  at  with 
axes  respectively  parallel  to  those  of  C^-xyz. 

(and  thus  also  parallel  to  those  of  A-xyz  ). 
a  Cartesian  coordinate  system  centered  at  J.  with 
axes  respectively  parallel  to  those  of  C^.  -xyz. 


Vector  Notation 


I -I  I 


a  6x1  vector  used  for  notational  brevity  (defined  on 
page  30). 

the  3x1  vector  of  Cartesian  coordinates  of  B 
relative  to  A-xyz. 

the  3x1  vector  of  Cartesian  coordinates  of  C^ 
relative  to  A-xyz. 


R„  , ,  -  the  3x1  vector  of  Cartesian  coordinates  of  C, 

Ci/Ji  1 

relative  to  iK-xyz. 

Rrl/1  -  the  3x1  vector  of  Cartesian  coordinates  of  C, 

Ci/Ji  1 

relative  to  J!-xyz. 

"R* •*  -  the  3x1  vector  of  Cartesian  coordinates  of  J, 

Ji  1 

relative  to  A-xyz. 

The  orientation  (for  a  typical  position)  of  some  of  the  points, 
coordinate  systems,  and  vectors  appears  in  Figure  1  (next  page). 


Matrix  Notation 


A. 

J 


‘‘■’V 


.I-II1(Tj) 


>i 


I, 


i/C, 


V/c. 


vi-n 

J 


a  6x3  matrix  used  for  notational  brevity. 

(T.  is  defined  below). 

J 

the  whole-body  inertia  tensor  (3x3)  about 
axes  of  A-xyz. 

the  inertia  tensor  of  component  i  about 
the  axes  of  A-xyz. 

the  inertia  tensor  of  component  i  about 
the  axes  of  C^-xyz. 

the  inertia  tensor  of  component  i  about 
the  axes  of  C ^  -xyz. 
a  6x6  matrix  that  is  a  function  of  the 
elements  of  T.  (explicitly  defined  in 

J 

Equation  (4-8),  page  30). 
a  3x3  matrix  which  contains  elements  result¬ 
ing  from  application  of  the  parallel  axis 
theorem,  (more  explicitly  defined  in 


Equations  (2-9)  and  (2-10),  page  10). 


Qj 


h 


I-II 


I-III 


-  a  12x6  matrix  used  for  notational  brevity 


1a 


I-II 


I-III 


.111 'ooo 


-  a  13x6  matrix  used  for  notational  brevity 


Analysis 

The  whole-body  mass  properties  (with  respect  to  A-xyz)  can  be  expressed 
in  terms  of  the  component  properties  (with  respect  to  A-xyz)  as  follows: 

n 


Mass 


mB  =  X  mi 


(2-1) 


_  11  f 

Center  of  Mass  Rr  =  X 


i=i 


Rr  .  m. 

Li  1 


Mr 


(2-2) 


n 

Inertia  Tensor  In  =  J  I,. 
-  8  1 


(2-3) 


But 


Rci  =  Rc1/Ji  +  Ro. 


(2-4) 


Also  Rp  /-i  ~  T . Rp  i  # . 

ci/Ji  i  ci/Ji 


where  the  transformation  matrix  T.  is  given  by: 


A 


cosxcosy  cosysinysinv-sinxcosM 

sinxcosy  sinxsinysinu+cosxcosu 

-siny  cosysiny 


(2-5) 


cosys i nycosy +s i nys i ny 
sinxsinycosy-cosxsiny 
COSyCOSy 


(2-6) 


IU 


In  the  above  matrix,  x>  Y>  and  y  are  rotation  angles  illustrated 
in  Figure  2.  In  the  field  of  flight  mechanics,  y ,  x>  and  y  correspond, 
respectively,  to  aircraft  bank,  heading,  and  flight  path  angles  (Ref.  13, 
page  44).  The  orientation  of  the  primed  axes  with  respect  to  the  unprimed 
can  be  described  in  terms  of  these  three  rotation  angles.  To  define  these 
rotations,  two  intermediate  coordinate  systems  are  introduced  whose  properties 
are  as  follows:  the  system  Ox^y^z^  is  obtained  from  the  unprimed  system  by 
means  of  a  rotation  x  around  the  z-axis;  the  system  Ox^y^z^  is  obtained 
from  Ox-|y^z^  by  means  of  a  rotation  y  around  the  y^-axis;  the  primed 
axes  system  is  obtained  from  0x2y2z2  by  means  of  a  rotation  y  around  the 
X2~axis. 

Substituting  Equation  (2-5)  into  Equation  (2-4)  and  then  substituting 
(2-4)  into  (2-2),  the  whole-body  center  of  mass  becomes: 

rb  =  I  <Ti 
“  1=1  1 

By  applying  the  parallel  axis  theorem,  it  is  possible  to  express  I.. 

in  terms  of  I . /r  as  follows: 

1/ui 

Ii  =  I./c_  +  P.  i=l,...,n  (2-8) 

where  P^.  is  a  3x3  matrix  whose  elements  are  defined  as  follows: 


Off  Diagonal 
Elements 

Pi(q»r) 


1-  Rr  (q)  •  Rr  (r) 

i  l  Lj  Li  J 


(2-10) 


where  q,  r  e  {1,2,3}  and  q  f  r. 

Now  the  inertia  tensor  can  be. rotated  by  a  similarity  transformation, 

1  -1 

i.e.,  I  =  TI  T  ,  where  T  is  the  transformation  matrix  from  the  primed 
to  the  unprimed  axes  systems  (Ref.  10,  page  370).  Since  our  trans¬ 
formation  is  orthogonal,  then  T-'1  =  TT,  where  the  T  superscript  indicates 
the  transpose  (Ref.  14,  page  207).  Thus  I  =  TI'JT. 


For  I./r  ,  then,  we  can  write: 
1/Li 


1  i/C-  =  Ti  li  VC.  Ti 


Substituting  (2-11)  into  (2-8),  we  obtain: 


h  =  Ti  ri  vc,Ti  +  pi 


Substituting  (2-12)  into  (2-3)  results  in: 


(2-11) 


(2-12) 


=  l  T.I. l/r  T.1  +  P. 

1=1  l.  11  /ci  1  1 


(2-13) 


Note  that  in  Equations  (2-7)  and  (2-13),  that  while  T. (i=1 , . . .  ,n) 

and  P. (i=l , . . . ,n)  are  functions  of  position,  Rr,,,  (i=1,...,n)  and 
1  Li/0i 

I  ' 

i / C ^  (i=1...,n)  are  not  (assuming  components  are  rigid  bodies).  Equations 
(2-7)  and  (2-13),  then,  allow  us  to  determine  the  whole-body  mass  properties 
for  any  position,  given  the  body  position  and'  component  mass  properties  with 
respect  to  component-fixed  axes  systems.  A  computer  program  capable  of 
accomplishing  this  (PROGRAM  WHOLE)  is  provided  in  Appendix  A. 


CHAPTER  3 


The  Hanavan  Model 

A  personalized  human  body  model  was  developed  by  Ernest  P.  Hanavan 
(Ref.  11),  who  obtained  approximations  for  the  Rr,,,  and  I . ,  /r  (i-1 . n) 

V  °i  1  /Li 

by  assuming  simple  geometric  shapes  for  the  body  components.  Figure  3 
illustrates  his  model  (from  Ref.  11,  page  8),  and  the  assumed  component  shapes 
are  indicated  below: 


1 

Head 

-  right  circular  ellipsoid  of 

revol ution . 

2 

Upper  Torso 

-  right  elliptical  cylinder. 

3 

Lower  Torso 

-  right  elliptical  cylinder. 

4,5 

Hands 

-  spheres. 

6,7 

Upper  Arms 

-  frustrums  of  right 

circular 

cones . 

8,9 

Forearms 

-  frustrums  of  right 

circular 

cones . 

10,11 

Upper  Legs 

-  frustrums  of  right 

circular 

cones. 

12,13 

Lower  Legs 

-  frustrums  of  right 

circular 

cones . 

14,15 

Feet 

-  frustrums  of  right 

circular 

cones . 

Using  25 

anthropometric 

body  measurements  of 

the  person  being  modeled 

(Ref.  11,  pages  9-10,  A-l  through  A -3 )  and  Barter's  regression  equations  for 
estimating  masses  of  components  (see  Chapter  6),  Hanavan  defined  his  simply 
shaped  components  and  provided  formulas  for  component  mass  properties 
(Ref.  11 ,  pages  13-24) . 

In  this  way,,  personalized  approximations  can  be  made  for  whole-body 
mass  properties  in  any  position.  Hanavan  compared  his  model 's  prediction 
with  experimental  data  (Ref.  15),  and  claims  his  center  of  mass 
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Figure  3.  Hanavan’s  Mode] 
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is  generally  within  0.7  inches  and  moments  of  inertia  are  generally  within 
10"  of  experimental  values  from  pendulum  tests  of  live  subjects.  For  prelim¬ 
inary  design  studies,  this  may  be  sufficiently  accurate.  Where  better  accuracy 
is  needed,  the  method  presented  by  the  present  study  (Chapter  4)  can  be  used. 

The  work  of  Hanavan  was  expanded  by  J.  Dubois  and  others  (Ref.  7) 
who  added  a  model  pressure  suit  consisting  of  "shells"  of  similar  shapes 
surrounding  each  component.  They  also  improved  Hanavan 's  original  model  by 
assuming  the  torso  hollow  (thus  taking  air-filled  lungs  into  account). 

Further  improvement  of  Hanavan's  work  was  made  by  Tieber  and  Lindemuth 
(Ref.  17).  Though  the  shape  of  Hanavan's  components  is  not  radically  altered, 
pivot  points  are  changed  (Ref.  17,  pages  16-21).  Also,  a  different  set  of 
regression  equations  is  used  for  estimating  component  masses  (discussed  in 
Chapter  6).  A  pressure  suit  model,  similar  to  that  of  Dubois,  is  also  pro¬ 
vided  (Ref.  17,  pages  23-34),  based  on  a  G-4c  suit  (used  in  the  Gemini 
project).  Tieber  and  Lindemuth  claimed  that  the  center  of  gravity  errors 
are  generally  reduced  to  less  than  0.3  inches,  but  the  moment  of  inertia 
errors  are  generally  as  high  as  15%  of  experimental  values. 


The  present  study  presents  a  more  general  method  for  obtaining 

K.,,  and  I.'/r  (i=l,...,n)  than  Hanavan’s  model  (Chapter  4).  But  this 
‘'i/Ji  1/Li 

rethod,  it  will  be  seen,  depends  on  obtaining  estimates  for  the  sums  of 


component  inertia  moments.  Hanavan's  model  provides  one  possible  method  of 
obtaining  these  estimates  (see  Chapter  7). 


This  more  general  method  will  also  depend  on  estimates  of  component 
r.isses.  Barter's  regression  equations  will  be  presented  in  Chapter  6  as  a 


sible  estimation  method. 


CHAPTER  4 


Determination  of  Component  Mass  Properties 
About  Component-Fixed  Axes 

Overview 

As  was  shown  in  Chapter  2,  whole-body  mass  properties  for  any  position 
can  be  found  if  component  mass  properties  about  component-fixed  axes  are  known. 
Hanavan's  model  of  the  previous  chapter  obtained  these  by  assuming  simple 
geometric  shapes  for  the  components.  This  chapter  provides  a  more  general 
method  which  assumes  nothing  about  the  component  shapes,  only  that  the 
components  are  rigid  bodies.  It  does  depend  on  the  assumption  that  the 
following  information  can  be  obtained: 

(1)  Experimental  data  for  the  whole-body  center  of  mass  and 

the  inertia  tensor  for  a  living  human  (discussed  later  in  this  chapter). 

(2)  An  estimate  for  the  mass  of  each  of  the  various  components 

of  the  body  (defined  by  the  researcher).  Chapter. 6  presents  two  methods, 
one  using  either  of  the  two  sets  of  regression  equations  referred  to 
in  Chapter  3,  the  other  using  experimental  data  for  casts  of  the 
particular  subject's  body  parts. 

(3)  An  estimate  for  the  sum  of  the  three  moments  of  inertia  for 
each  body  component.  Chapter  7  provides  two  methods,  one  using  Hanavan's 
simply  shaped  components,  the  other  using  experimental  data  for  casts 

of  the  subject's  body  parts. 

The  basis  for  the  solution  for  component  mass  properties  is  the  formation 
c/  what  will  be  called  mass  property  difference  equations.  These  are  formed 
t/  subtracting  the  equations  for  whole-body  mass  properties  (Equations  (2-7) 
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and  (2-13))  for  selected  body  positions.  Using  the  center  of  mass  difference 
equation  and  estimates  for  component  masses,  component  centers  of  mass  can 
be  determined.  Using  the  inertia  tensor  difference  equation,  estimates  for 
component  masses,  and  estimates  for  sums  of  component  moments  of  inertia, 
component  inertia  tensors  can  be  determined.  It  will  be  shown  in  Chapter  5 
that  estimates  for  sums  of  component  moments  of  inertia  provide  a  needed 
additional  independent  equation  relating  elements  of  component  inertia  tensors. 


Component  Centers  of  Mass 

For  two  different  body  positions,  say  I  and  II,  identical  except  for  the 
position  of  one  component,  say  number  j,  Equation  (2-7)  becomes: 

f  n 


R 


I 


B 


l  (Ti  !R  +  Rj  )  ■  rni 
i  =  l  1  c!  /J.  Ji  1 


/m 


B 


R  11 
KB 


f  <h%c,  /0.  +R”)  -J 

(i  =  i  J  i  '  i  J 


/m 


B 


Splitting  up  the  summation  and  subtracting  these  two  equations: 

r  n  T  n  t  T  n  T  n 


«bt  -V1 


Equivalently: 


I-II/k- 


‘  A  ’  (RB}  =  E 
B  B  i=l 


Rr  1  , ,  m .  + 
Ci/Ji  l 


a1-11!?,  ) 

Ji  j 


(4-1) 


But  I  and  II  are  identical  positions  except  for  the  position  of  component 


Thus : 


10 


Ti  =  1=1.  ...»  n  (i  f  j) 

RJ.  =  RJ!  i=l,  ....  n 

Hence  Equation  (4-1)  becomes: 

mB  =  [4l'n<y]  rc,/j.  mj 

j  j 

Equivalently: 

Al-II(Jj)  RC1/J<  =  )  aI-II(Rb)  (4-2) 

J  J 

In  the  above  equation  Rg  and  Rg1  can  be  determined  by  experiment 

(see  page  ),  as  can  mg  .  The  quantity  m-  can  be  estimated  (see 

page  ).  The  matrices  T^  and  T^1  represent  arbitrarily  selected,  known 

J  J 

positions.  Thus  Equation  (4-2)  represents  a  system  of  three  linear  equations 

—  J  I 

in  three  unknowns,  the  elements  of  Rri/1  .  Clearly,  if  A  (T.)  is 

nonsingular,  Rr , , ,  can  be  solved  for  directly  in  Equation  (4-2).  However, 
o./J. 

J  J  J  T  J 

the  following  analysis  proves  that  A  (T.)  is  singular  for  all  choices  of 

J 

positions  I  and  II. 

The  orientation  of  the  component-fixed  axes  Cj-xyz  1S  arbitrary.  Thus, 
without  loss  of  generality,  C’.-xyz  can  be  oriented  such  that  x*  =  0°, 

J 

y 1  =  0°,  and  y1  =  0°  for  position  I,  for  which: 

1  0  0 
0  1  0 
0  0  1 


For  any  other  position,  II: 
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IH  = 


_  II0  II 
COSX  COSY 


sinx^cosy11 


-siny11 


cosx^siny^siny** 

-sinx^cosy^ 

sinx^siny^siny1* 

A  II  II 
+cosx  cosy 

cosy^siny* 1 


„nc  II.  II  II 
COSx  Siny  cosy 

.  II  .  II 
+sinx  siny 

sinx^siny^cosy^ 

II  .  II 
-cosx  siny 

II  II 
cosy  cosy 


Hence: 


.1-11 


'V 


,  II  II 
1-COSx  cosy 


.  II  II 
-siny  cosy 


sin,11 


-cosx^siny^siny*1 

+sinxIIcosy11 

,  -  II  -  II  -  II 

1-smx  siny  siny 

-COSy^COSy  ^ 

-cosy^siny^ 


II  .  II  II 
-cosx  siny  cosy 

.  II  .  II 
-sinx  siny 

-sinx^siny^cosy^ 

II  .  II 
cosx  siny 

1  -cosy^cosy^ 


The  determinent  of  the  above  matrix  turns  out  to  be  zero.  This 


U  seen  by  expanding  the  determinent  and  substituting  appropriate  trigono- 
'irtric  identities.  Consequently,  for  any  selection  of  positions  I  and  II, 
matrix  A  **(T.)  is  singular. 

J 

Hence  Equation  (4-2)  represents  a  dependent  system  of  equations 

"'  ■  *  1^»  133),  and  therefore  does  not  uniquely  determine  the 

«!<?«•:*  t$  of  .  However,  if  an  arbitrary  third  position  (III)  is 

J  0 

!**;  f,cod,  the  following  difference  equation  results,  identical  in  form  to 


•r.Scn  (4-2): 


Rc:/j.  =  <W 
0  J 


H-!n(RB) 


(4-3) 


combination  of  Equations  (4-2)  and  (4-3)  constitute  a  system  of 
*  Motions  in  three  unknowns.  For  judicious  selection  of  positions  I, 
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II,  and  III,  it  turns  out  that  this  system  can  contain  three  independent 

equations,  thus  determining  a  unique  solution  for  Rpl/1  .  However,  it  may 

J7  J 

be  inconvenient  to  isolate  three  independent  equations.  An  alternative 
solution  of  the  system  which  does  not  require  this  isolation  does  exist,  and 
its  derivation  follows.  • 

Combining  Equations  (4-2)  and  (4-3)  into  one  equation: 


AI-II(Tj) 

A1'111^) 


J 


aMI(Rb) 


aI-hi/o  a 
a  (M 


Defining  the  matrix  A.  = 

J 


equation  by  A^.: 


and  multiplying  the  above 


AjAjRCl/J.  =  Aj 

J  J 


(4-4) 


But  for  any  matrix  A-,  A'-A.  and  A.  have  the  same  rank,  or  the 

J  J  J  J 

number  of  independent  rows  (or  columns)  (Ref.  14,  page  139).  Thus 
if  A,  has  rank  3,  then  so  has  A^A-.  But  A^A-  is  3x3;  hence,  aTa. 
rust  have  an  inverse  if  it  has  rank  3.  Therefore,  if  A.  is  of  rank  3, 

J 

then  Equation  (4-4)  can  be  solved  directly  for  R~ ,  , ,  as  follows: 

Cj/Jj 


Rc:/J,  -  (m  /m.)  (aTa.)"^aT 

J  J  D  J  J  J  J 


A’-n(RB) 


(4-5) 


Whether  or  not  A.  has  rank  3  depends  upon  the  selection  of  positions 
J 

*.  II,  and  III.  Two  specific  examples  follow,  the  first  of  which  is  an  A. 

3 

s'  rank  3,  the  second  of  which  is  an  A.  of  rank  2,  thus  substantiating  the 

J 

<Jrlier  statement  that  the  selection  of  body  positions  must  be  judicious. 
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Example  1 

Let  positions  I,  II,  III  have  rotation  angles  as  indicated  below. 


POSITION 

X . 

Y 

u 

I 

0 

0 

0 

II 

0 

sin"1 (3/5) 

0 

III 

0 

0 

sin"1 (3/5) 

For  this  selection  of  positions: 


Thus: 


Y 

0 

o' 

"  4/5 

0 

3/5' 

"l 

0 

0  ‘ 

0 

1 

0 

Tj1* 

0 

1 

0 

Tjn  “ 

0 

4/5 

-3/5 

0 

0 

1_ 

-3/5 

0 

4/5_ 

0 

3/5 

4/5_ 

V 


^-"(V 

?:Tir(V 


1/5 

0 

-3/5 

0 

0 

0 

3/5 

0 

1/5 

0 

0 

0 

0 

1/5 

3/5 

0 

-3/5 

!/5_ 

For  the  above  matrix  A.: 

J 


0 


0.4  0 

0  0.4  0 

0  0  0.8 


is  clearly  of  full  rank.  Thus  (aTa.)"1  exists  and  is  in  fact  as 

J  J 


f  ,  M  , 


(ATA.)-1 
v  J  Y 


2.5 

0 

0 


0  0 
2.5  0 

0  1 .25 


Therefore,  for  this  selection  of  body  positions  Equation  (4-7)  will 
provide  a  unique  solution  for  the  elements  of  •  In  Example  2,  which 

follows,  this  will  not  be  the  case. 


Example  2 

For  this  example  let  positions  I,  II,  III  have  rotation  angles  as 
indicated  below. 


POSITION 

X 

Y 

y 

I 

0 

0 

0 

II 

0 

sin"1 (3/5) 

0 

III 

0 

sin"1 (4/5) 

0 

For  this  selection  of  positions: 


'1 

0 

0' 

'  4/5 

0 

3/5' 

‘  3/5 

0 

4/5" 

0 

1 

0 

1’  = 

0 

1 

0 

■ 

0 

1 

0 

0 

0 

l 

_-3/5 

0 

4/5_ 

_-4/5 

0 

3/5_ 

1/5 

0 

3/5 

0 

0 

0 

-3/5 

0 

1/5 

2/5 

0 

-4/5 

0 

0 

0 

_  4/5 

0 

2/5 

fc.r  the  above  matrix  A: 


*1*. . 

1.2  0  0  “ 

0  0  0 

0  J 

0  0  1 .2 

which  is  clearly  of  rank  2.  Therefore,  A(A.  will  not  have  an  inverse. 

J  J 

The  matrix  A.  provides  only  two  independent  equations  and  thus  does  not 
.  3 

determine  a  unique  value  for  the  elements  of  Rrl/,  . 

Cj/Jj 

An  apparent  (though  not  proved)  general  rule  for  the  selection  of  body 
positions  to  insure  a  unique  value  for  is  to  choose  positions  such  that 


III 


J  3 


R-1,.  >  Rr^/i  >  and  R*1/.  are  not  all  three  coplanar.  Thus,  in 
YJi  Ci/Ji  Ci/Jj 

Example  1,  where  these  three  vectors  did  not  lie  in  the  same  plane,  a  unique 


value  for  Rr ,  , ,  was  determined.  But  for  Example  2,  where  the  vectors 
Y  3 

were  coplanar,  a  unique  value  was  not  determined. 

In  the  analysis  of  this  section,  it  was  tacitly  assumed  that  for 
component  j,  it  is  possible  to  define  positions  I,  II,  and  III  such  that 
only  j  is  different.  But  this  is  not  possible  for  some  components,  for 
example  if  the  thigh  position  is  different,  then  so  must  the  calf  and  foot 
positions  be  different.  More  generally,  Equation  (4-5)  is  applicable  directly 
vA y  to  those  components  connected  to  only  one  other  component,  hereafter 
Silled  external  components.  Those  components  which  are  connected  to  two 
r.ber  components  shall  be  called  internal  components.  Note  that  the  central 
component  from  which  the  limbs  emanate  fit  into  neither  of  these  categories 
s*  components,  as  more  than  two  components  are  connected  to  it.  For  conven- 
this  central  component  will  be  numbered  n  .  Typically,  of  course, 
f'’ ;  central  component  is  simply  the  torso,  but  the  researcher  has  the  option 
U  fr elude  the  head  and  neck  as  part  of  the  central  component  if  they  are 
* < ••  :  to  be  immobile.  This  assumption  may  be  quite  reasonable  for  a 

V re- suited  astronaut,  with  restricted  head  mobility. 

'  method  for  determining  the  center  of  mass  for  any  internal  component 
s  1  :"</  one  component  external  to  it  is  now  presented.  Using  the  method 


1 ' * presented,  the  center  of  mass  of  the  external  component  can  be 


determined.  For  notations!  convenience  let  this  external  joint  be  numbered 
j+1  and  the  internal  component  whose  center  of  mass  we  seek  be  numbered  j. 
Then  difference  equations  similar  to  Equations  (4-2)  and  (4-3)  can  be  formed 
but  will  have  an  extra  term  (because  the  contribution  of  component  j+1 
does  not  subtract  out)  and  will  be  as  follows: 

_iI'II(Tj)]  V./J. c  . 

-  P'n<wvj+1  +  iI'II(V  " 


a1 -i 1 1 (t  .  n  rci/j 

L  j  J  Cj/j j 


mB/mj  ^-in(RB) 


mj+-]/m-  A  (Tj+-|)R0'  /j  +  A  (Rj  ) 
3 1  L  3  '  Lj+l/dj+l  Jj+1 


In  the  above  equations,  the  center  of  mass  of  the  external  component, 

,,  must  be  predetermined.  The  T.,  T.,,,  and  R,  are 

Cj+1/Jj+1  J  J+1  Jj+1 

functions  of  the  body  positions  I,  II,  and  III.  Combining  the  above 

equations,  an  equation  analagous  to  Equation  (4-5)  is  obtained  and  is  as 

follows: 


f  T  r1  T  r  &MI<0 
Rf/i  “  A  .A  •  A.  mR/m.  -r-Trr-- — 
Cj/Jj  J  B  ai"iii(rb) 


.  ,  4  (Ti+, )  _  4  <RJ  ) 

-  m.  ,/m,  J - -—  Rr,  71  +  r _ _i±l _ _ 

J 1  j  n  ai-hi(t  cj+i/Jj+i  i_ni,_  s 
LA  lij+r  a1  ) 

L  L  dj+i 


A  similar  analysis  can  be  carried  out  to  find  the  center  of  mass  for 
each  internal  component  with  more  than  one  external  component.  The  following 


equation  results: 


—  f  t  I  T 

C'./Jj  to  jJ  J 


l  |V'n(RB) 


i  external  A  vT - )  I"UI(r  ) 

to  j  1  L  o. 


The  above  equation  can  actually  be  considered  generally  applicable  to 
any  component  except  the  central  component  if  the  summations  are  ignored 
for  external  components  (for  which  there  are  no  i  external  to  j  ).  Assum¬ 
ing  Rr,/-,  (i=l .  ,n-l )  have  already  been  determined,  the  center  of  mass 
Li/Ui 

of  the  central  component  is  given  by  a  simple  manipulation  of  Equation  (2-7) 
for  any  position: 


Lr ,,  =  T  (m„/m  )  FL 
C^/J  n  B  n  B 


n 


1  -  R- 


Jn  in  the  above  equation  can  be  any  conveniently  defined  fixed  point 
cr,  the  central  component. 

Using  the  equations  developed,  then,  all  component  centers  of  mass  can 
found. 


“ "■|,T;'r'Ont  Inertia  Tensors 

for  two  different  body  positions  I  and  II,  identical  except  for 
tv  position  of  one  external  component  j,  Equation  (2-13)  becomes: 
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I\T  .  A) 


l  Ti  xi/c <T> +  pi 


i 


ii 

B 


i=n  1  i/ci 


1R1  ii/c/T!I)T  +  p!n 


i=l v  •  '  i 

Subtracting  these  two  equations: 


n  / 


AI_II(i  )  =  i1  -i11  =  y  Mi.  ( t  1 ) T-j 1 1 1  *  (t1  1  )T+pI-pI I] 

A  lV  B  .  i=1l  i  i/C1*  i;  1  i/C^i  '  i  i 


-1!  i  • 

i  ■ 

But  for  all  components  except  j, 


T>  -  T” 

9 

i=l , . . .  ,n 

(i7j) 

P;  =  p” 

9 

i=l . .  ,n 

(i7j) 

Substituting,  the 

following  equation 

is  obtained 

(tb  ) 


(T]1)1 


Equivalently: 

*MI<Ib>  -  ■ T]  ijVcjtTj)7  -  T"  Ij^c/rj1)1  '4-6> 

In  this  last  equation,  Ig  and  1^  are  obtainable  from  whole-body 
rental  data  (see  Chapter  6).  As  can  be  seen  in  Equations  (2-9)  and 


M3),  P.  is  a  function  of 

J 

V 

T-,  PVI/1  ,  and  R 

J  Cj/Jj  Ji 

The  choice  of 

<:  f>  positions  determines  T  . 

J 

and 

Rj  /A’  The  QuaHti ty  itu 

3 

is  obtainable 

< •  i-t'ore  from  the  methods  presented  in  Chapter  6.  The  quantity  is 


i”:  ‘-i •'  s!>le  from  the  method  presented  in  the  previous  section.  Thus  P.  can 

3 

■  •  '  for  any  positions  I  and  II,  and  Equation  (4-6)  becomes  a  matrix 

in  one  unknown  matrix,  I-'/r  .  However,  there  is  no  obvious  way  to 

J/  J 

•  ‘  > *  (i- 0)  directly  for  I^y^  .  Hence  it  becomes  necessary  to  express  the 
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matrices  in  terms  of  elements,  carry  out  the  operations,  and  equate  elements. 

In  this  way  Equation  (4-6)  is  equivalent  to  the  following  system  of  nine 

linear  equations  in  nine  unknowns  (the  elements  of  !j/C,  ): 

J 


Main  Diagonal  Elements 


3  3  r  i 

AI_II(lB(pJp))-AI'II(Pi(p,p))  =  l  l  AI-II(T.(q,k)T.(r ,£,))  I  '  (k,<0 

B  J  <1=1  k=lL  J  J  J 


where  p,q,r  e  (1,2,3)  (p/q^r) 


Off  Diagonal  Elements 


3  3 

‘  '(IB(q.'-))-4I'II(P1.(q.r))  -  I  I  [J-n(T.(q,k)T1.(r,t))]I1)c  (k.t) 

J  £=1  k=1  J  J 

where  q,r  e  (1 ,2,3}  ,  q^r 

Because  of  the  symmetry  of  the  inertia  tensor,  the  above  system  can  be 
nduced  to  a  system  of  six  equations  in  six  unknowns  (the  six  independent 
fluents  of  I . )r  ).  The  off  diagonal  equations  given  above  actually  provide 

-  j 

:»•!/  three  rather  than  six  equations.  For  example  if  q=l ,  r=2,  then  above 

tuition  becomes: 


3  3 

AI-II(In(l ,2))-AI_II(P-(l  ,2))  =  l  l  [a(T.(1 ,k)T. (2  ,$,))]  I,'  (k,i 

J  £=1  k=1  J  J  J/tj 

But  for  q=2,  r=l  : 


3  3 

4'n<IB(2,l))-4I'II(P.(2,l))  =  M  [4(T,(2,k)l\(l,t))]  I.'  (k.t) 

J  £_■]  |^_  "J  J  J 

tv*  Ip  >  iP  >  pL  P^ 1 ,  and  I  '  are  symmetric  matrices,  or: 

A  o  J  .J  J/Uj 


ij  (2,1)  =  Ig  0,2) 

Ig1  (2,1 )  =  Ig1  (1,2) 

Pj (2 ,1 )  =  Pj ( 1 ,2) 

P j 1 ( 2 ,1 )  =  P j 1 ( 1 ,2) 

■I0/Cj(k,t)  =  Ij/Cj(^k) 

Substituting  these  into  the  earlier  equation  for  q=2,  r=l,  reversing 
the  order  of  summation,  and  commuting  the  multiplications  within  the  summa¬ 
tion: 


2'!I(IbO,2))-4I‘II(P1(1,2))  -  l  l 

B  J  k=l  Z=1 

But  k  and  i  are  dummy  variables.  Hence  the  equation  above  for 

;r2,  r=l  is  seen  to  be  identical  to  the  earlier  presented  equation  for 

■V 1,  r=2.  Similarly,  identical  equations  result  for  q=l ,  r=3,  and  for 

1"3,  r=l.  The  same  is  true  for  the  cases  q=2,  r=3  and  q=3,  r=2. 


AI"II(Tj(2,k)Tj(l,£)) 


In  this  way  the  original  system  of  nine  equations  reduces  to  the  follow- 
linear  system  of  six  equations  in  six  unknowns  (the  independent  elements 

i »  j  >  )  • 

jV 

Kain  Diagonal  Elements 


i*’;i(Ip(p,p))-AI’II(Pi(p,p))  =  l 


k=l 


AI-II(Tj(p,k)^) 


^VC^P’P) 


3  3  r 

+2  l  l 

£=1  k=l 
k>£ 


a(Tj (q ,k)Tj (r ,k) ) 


JIj/c.(£,k) 

J 


where  p,q,r  {1,2,3}, 


p/q^r  q<r 


Off  Diagonal  Elements 


(IB  (q,r))-4I-II(Pi(q,r))  =  'l  4I'II(Tj(q,k)Tj(r,k))  I^c  (p,p) 

k~l  J  J 


+  l  l  [4I'II(Ti(q,k)Ti(r,«) 
£=1  k=l  3  J 

k>£ 


+  Tj(q,i)Tj(r*»k)Oijyc>u,k) 

where  q,r  e  {1 ,2,3} ,  q<r 

In  matrix  form  this  linear  system  can  be  written  as  follows: 
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l 

.if 

>;■ 


% 


where 


and 


b^"11  = 


A1"”  (Ig  0,1)) 

A1"11  (IB  (2,2)) 
A1"11  (Ig  (3,3)) 
A1"11  (Ig  0,2)) 
A1'11  (Ig  0,3)) 
A1"11  (lB  (2,3)) 


A1-11  (Pj  (1,1)) 
AI"II(Pj(2,2)) 

AI-II(Pj  (3*3.) ) 

AI_II(Pj(l  ,2)) 

A1-11 (Pj (1  ,3) ) 
AI“II(Pj(2,3)) 


l-U 

J 


r-n 


AtTjd.i)2) 

4(Tj(1.2)2) 

A(rj0.3)2) 

2(»TJ(1.1)Tj(1.2)) 

2(aTj(1 ,1 )Tj(1 ,3) ) 

2(ATJ(1.2>TJ(1,3)) 

i(Tj(2.1)2) 

A{Tj{2.Z)2) 

A(Tj(2.3)2) 

2(ATJ(2.1)TJ(2,2)) 

2<4Tj(2.1)Tj(2.3)) 

2(ATJ(2,2)Tj(2.3)) 

AtTjO.n2) 

4(Tj(3,2)2) 

a(Tj(3,3)2) 

2(aTj(3 ,1 )Tj(3.2) 

2(«TJ(3.t)Tj(3,3)) 

2(aTj(3,2)Tj(3.3)) 

a(Tj(2.1)Tj(3,I)) 

a(Tj(2.2)Tj(3.2)) 

A (Tj (2  »3)Tj (3  *3) ) 

ri(TJ(2.2)TJ(3.i»1 

[♦«(Tj(2.))Tj(3,2))J 

[  4{Tj(2,3)Tj(3,1))1 
1  L*4(Tj(2.1)Tj(3.3))j 

'  4(TJ(3.3)Tj{3.2))' 
♦i(Tj(2.2)T1(3,3)), 

a(Tj(1 .l)Tj(3,l)) 

a(Tj(1 ,3)Tj(3.3)) 

r  4(TjO. 2)Tj(3. 0)1 
U(TjO.  DTjO.  2))J 

1  [  4(Tj(1.3)Tj(3.1))‘ 
[+4(Tj(l  ,3 )Tj(-3»3) }j 

1  [  A(Tj(1.3)Tj{3,2))' 
[,4(Tj(1.2)Tj)3.3)). 

L 

A(Tj(1.2)Tj(2,2)) 

4(Tj(1.3)Tj(2,3)) 

[  s(Tj{1.2)TJ(2.1))] 
L,A!TJU,mj(2.2)). 

r  A(TJn,3)TJ(21i))| 

I  L+AfTjO  .3  )Tj(2.3)), 

'  4(TJ(l,3)Tj(2.2))‘ 
+a(Tj(1  .2)Tj(2.3)J. 

Clearly,  if  the  6x6  matrix  V. 

V 


I-II 


(4-8) 

has  an  inverse,  then  the  vector 


ekc.ents  of  I^g  can  be  solved  for  directly  in  Equation  (4-7). 


\  it  turns  out  that  V- 

J 


I-II 


is  singular  for  any  selection  of  positions 


5;  rigorously  proved  in  Chapter  5).  In  fact,  for  every  one  of  numerous 

is  only  4,  even  when  position  II  's 
’  .  arc  all  different  than  position  I's. 

*  third  position  (III),  identical  to  I  and  II  except  for  the  position 


:  c/anples,  the  rank  of 

J 
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of  component  j,  is  introduced,  then  an  equation  like  (4-7)  can  be  written, 
as  follows :  ' 


V1'111 


’j/c'1-1) 

J 

hk™ 

luc™ 

^•/C^1  ’2) 

lik"-3) 

Tj/C^2,3) 


=  b 


I- 1 II 


(4-9) 


.  ,.1-m 

where  V. 

J 


and  bf 111  have  the  same  form  as  the  previously  given 

J 


equations  for  V?-**  and  b^  respectively. 

J  J 

Combining  Equations  (4-7)  and  (4-9): 


v  I-H  1 
J 

3 

V  I~111 

(2,2) 

= 

-  I-III 

b 

J 

1 

j 

(3,3) 

' 

0,2) 

0,3) 

(2,3)  . 

For  notational  brevity,  the  12  x  6  matrix 


I-II 

J _ 

i-iii 


(4-10) 


shall  be  denoted 

the  6x6  matrix 


ho  symbol  Q..  Now,  if  Q.  be  of  rank  6,  then  so  wil 

U  J 

.  In  this  case  (qTq.)""'  will  exist,  and  Equation  (4-10)  can  be  solved 

*3  J 

ell y  for  the  vector  of  the  elements  of  I  .'/r  .  However  it  turns  out  that 

J/ j 
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the  rank  of  Q.  is  less  than  6,  the  actual  rank  depending  on  the  selection 

J 

of  positions  I,  II,  and  III.  That  the  rank  of  Q.  is  less  than  six  is 

J 

proved  in  Chapter  5,  where  it  is  additionally  proved  that  no  matter  how  many 

positions  are  incorporated  into  Equation  (4-10),  the  Q.  matrix  will  still 

J 

be  of  rank  less  than  6. 

For  three  positions,  it  is  possible  to  find  specific  cases  for  which 
the  rank  of  Q.  is  5.  Thus  there  is  hope  for  a  unique  solution  if  one 

\J 

additional  equation  in  the  elements  of  I.'/r  is  introduced.  An  approxi- 

J/  J 

nation  for  the  sum  of  the  moments  of  inertia  of  component  j  provides  this 
equation: 


bvc^-^'j/c/2’2' 


+  'j/C  <3’3)  =  sj 
J 


(4-11) 


The  constant  s.  can  be  approximated  using  the  methods  presented  in 

J 

Chapter  7- 

Incorporating  Equation  (4-11)  into  (4-10): 


1* 

i 

r  ij)c.(i.D  j 

[ 

1 

* 

l,',c  (2,2) 

1  J/  0 

’d/c/3-3’  1 

= 

7—  I-III 
b  . 

1  J 

i 

%  ' 

L 

■I, 

! 

Ij/cd1,3) 

W:  ■ 

! . 

Jr  it  4 

bvcj12’3) 

(4-12) 
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i-n 

_CL___ . 

V1:111 
_J . . 

1110  0  0 

*  *T  *  _■] 

Now  if  Q  is  of  rank  6,  then  (Q  Q  )  exists  and  Equation  (4-12) 
can  be  solved  for  the  elements  of  Ij/Q,  as  follows: 

b/c/1-1' 
b/c/2-2> 

’j/ci3,3) 

Ij/Cj(1’2) 

1  j/Cjt1 ’3> 

b/c/2’3' 

An  example  of  body  positions  I,  II,  and  III  for  which  (4-13) 

*k-lds  a  unique  solution  for  the  elements  of  I  '  r  follows. 

J/Lj 

Example  3 

Let  positions  I,  II,  and  III  be  defined  the  same  as  they  were  for 
1.  For  these  positions. 
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j 

"  0.36 

0 

-0.36 

0 

-0.96 

0 

1 

i 

I 

0 

0 

0 

0 

0 

0 

-0.36 

0 

0.36 

0 

0.96 

0 

j 

0 

0 

0 

0.6 

0 

0.2 

1 

* 

0.48 

0 

-0.48 

0 

0.72 

0 

* 

f 

-t. 

0 

0 

0 

0.2 

0 

-0.6 

* 

qJ  ■ 

0 

0 

0 

0 

0 

0 

i 

0 

0.36 

-0.36 

0 

0 

0.96 

§ 

f. 

j 

0 

-0.36 

0.36 

0 

0 

-0.96 

ai' 

■c 

0 

-0.48 

0.48 

0 

0 

0.72 

% 

I 

r  ■ 

0 

0 

0 

o 

1 

0.2 

0 

$ 

| 

0 

0 

0 

0.2 

0.6 

0 

1 

1 

1 

0 

0 

0 

J 

m 

*sr  which 

"V 

.5$ 

;?y 

1.4896 

1 

0.5104 

0 

-0.3456 

0 

1 

1 

1.4896 

0.5104 

0 

0 

0.3456 

I 

*T  * 

Q,  Q,-= 

0.5104 

0.5104 

1.9792 

0 

0.3456 

-0.3456 

| 

J  J 

0 

0 

0 

1.52 

0 

0 

;•  > 

A  U 

-0.3456 

0 

0.3456 

0 

2.7616 

0 

: 

'  J'fc;1. 

0 

0.3456 

-0.3456 

0 

0 

2.7616 

M. 

■■Sv 

is,  it  turns  out,  of  full 

rank. 

and  thus 

*T  * 

(Qj  V 

exists. 

Thus 

$■■■  f iu rt ion  (4-13)  will  determine  a  unique  set  of  values  for  the  elements  of 


ft  i'.jaMon  (4-13)  is  not  applicable  to  internal  components  directly.  An 

f)  v.‘t'  *  similar  to  that  for  the  center  of  mass  determination  of  the  previous 

>;i_  ■ 

■«  :';<*  results  in  the  following  equation  for  the  inertia  tensor  elements  of 

'-jSsC  ■ '  ■ 

W’ 
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’Ms  last  equation  was  obtained  by  simple  manipulation  of  Equation  (2-13). 

The  last  two  sections,  then,  haye  provided  a  means  for  finding  the 
:-?r.tcr  of  mass  location  and  inertia  tensor  for  each  body  component  about 
exponent-fixed  axes.  A  computer  program  capable  of  using  this  means 
M'vGRAM  PART)  is  provided  in  Appendix  B.  Once  the  component  properties 
found,  Equations  (2-7)  and  (2-13)  can  be  used  to  predict  whole-body  mass 
fr.p?rties  for  any  body  position  of  a  particular  subject. 

f -l rental  Determination  of  Whole-Body  Mass  Properties 

At  the  beginning  of  this  chapter,  it  was  seen  that  experimental  data  for 
*<V*V-body  mass  properties  for  selected  positions  was  needed  in  order  to  apply 
‘V?  nethod  developed  by  this  chapter  for  prediction  of  component  mass  prop- 
This  section  briefly  describes  how  this  whole-body  mass  property 
right  be  obtained. 

Vhile  several  researchers  have  experimentally  determined  centers  of 
for  living  subjects,  only  one  determination  of  whole-body  center  of 
;rs '"*/  and  moments  of  inertia  has  been  made  (Santschi,  DuBois,  and  Omoto). 
•Via',*,  whole-body  products  of  inertia  have  not  been  experimentally  deter- 

However,  Fowler  provides  a  means  for  carrying  out  such  a  determination 
Santschi ,  DuBois,  and  Omoto  (Ref.  15)  experimentally  determined  centers 
«f' fXdC.y  and  moments  of  inertia  about  three  axes  for  66  living  subjects  in 
VV  -editions .  This  was  done  by  finding  oscillation  frequencies  of  a 
VH'Amf  pendulum  suspended  alternately  from  two  parallel  axes.  Mean  values 
•■*#  t't* tird  deviations  were  determined  for  the  subjects  studied.  Results 
:;>tr  '  V  at  useful  predictions  of  moments  of  inertia  for  similar  subjects 

-  V  • : V-  from  height  and  weight  alone. 

-v?  theoretical  accuracy  of  the  pendulum  used  was  about  0.5  percent 


for  the  center  of  gravity  and  between  2  and  .  percent  for  the  moments  of 
inertia.  The  length  of  the  pendulum  was  long  and  the  pendulum  oscillations 
stall  (1°)  to  minimize  the  effect  of  shifting  body  fluids. 

Fowler  (Ref.  9)  developed  a  theoretical  method  for  experimentally 
determining  spacecraft  mass  distribution  while  in  space.  By  using  a  series 
:f  three  thrusting  maneuvers,  the  center  of  mass  and  inertia  tensor  of  the 
craft  can  be  determined,  assuming  the  spacecraft  mass  is  predetermined  from 
i  linear  acceleration  and  that  angular  rates  can  be  measured.  This  method 
applicable  to  a  human  subject.  If  the  experiment  is  done  on  earth,  then 
v  :vn  torques  can  replace  the  known  thrusts.  Though  there  is  currently  no 
twice  available  with  which  to  spin  a  human  about  three  perpendicular  axes 
,>Mch  Fowler's  method  requires),  one  can  be  developed. 

In  order  to  minimize  the  effect  of  redistribution  of  body  fluids, 
im,wlar  rates  should  be  kept  as  small  as  is  accurately  measurable. 

If  it  is  desired  to  include  outer  coverings  in  the  values  for  mass 
,  then  the  subject  needs  only  to  wear  said  coverings  during  the 
-{rental  tests  (for  either  of  the  above  described  methods). 

I  .*■*;! ",  Ration  of  Outer  Garments  or  Suit 

'cr  some  application,  e.g.  a  maneuvering  astronaut,  the  subject  may 
r  outer  covering  that  moves  with  the  body  parts.  The  contribution  of 
v-w  t’.'tr  covering  to  the  system  mass  properties  may  not  be  negligible. 

: attribution  of  these  outer  garments  or  suit  can  be  accounted  for  by  using 
c \ ons  already  developed  if  we  consider  the  outer  covering  of  component 
*•  ;-t't  of  j.  No  changes  occur  in  the  equations  already  given  except 
f  1  and  s.  will  include  both  the  body  part  and  outer  covering  contri- 

■  *  *J 

m  It  is  shown  in  Chapters  6  and  7  how  m.  and  s.  estimates 

J  J 


'eluding  the  outer  covering  can  be  made.  This  method  will  be  valid  if  the 
'blowing  two  assumptions  are  made: 

(1)  The  components  of  the  outer  covering  corresponding  to  the 
body  components  are  themselves  rigid  bodies. 

(2)  There  is  negligible  relative  motion  between  each  component 
and  its  outer  covering. 

The  basis  for  the  solution  for  component  mass  properties  is  the  forma- 
!  ;r,  of  what  will  be  called  mass  property  difference  equations.  These  are 
V-ned  by  subtracting  the  equations  for  whole-body  mass  properties  (Equations 
M)  and  (2-13))  for  selected  body  positions.  Using  the  center  of  mass 
inference  equation  and  estimates  for  component  masses,  component  centers 
r'  asis  can  be  determined.  Using  the  inertia  tensor  difference  equation, 
for  component  masses,  and  estimates  for  sums  of  component  moments 
tM-'crtia,  component  inertia  tensors  can  be  determined.  It  will  be  shown 
"V'/.er  5  that  estimates  for  sums  of  component  moments  of  inertia  provide 
*  '.'  j  additional  independent  equation  relating  elements  of  component 

'  "WlM  tensors. 


CHAPTER  5 


Explanation  of  Infinity  of  Solutions 
for  Component  Inertia  Tensors 


In  Chapter  4  it  was  claimed  that  the  component  inertia  tensor  could 
ttt  be  uniquely  determined  without  the  introduction  of  an  additional  equation 
Siting  the  elements  of  the  inertia  tensor  (the  sum  of  component  moments  of 
’wrtia).  The  purpose  of  this  chapter  is  to  prove  this  claim. 

Equation  (4-8)  can  be  expanded  to  include  any  number  of  positions  as 
'.Allows: 


'ijf 

| 

A 

% 

f 

i'i 

■" 


& 


f 

V? 


w 

‘.ft  ‘ 

•k- 

m- 


rft, ' 

St 


1 

» — « 

M 

1 

I — t  mr-} 
> 

( _ _ 

Vi."111 

J 

Vl-N 

„  J  __ 

- 

V/c/1’1* 

(2,2) 

(3.3) 

(1,2) 

(1.3) 


r  ri-n  ] 
bj 

b1"111 

0 

- 

kN  J 

(5-1) 


5Qj 

N  is  the  number  of  body  positions. 

In  Chapter  4  it  was  claimed  that  the  inclusion  of  any  number  of  body 
y?.  will  not  provide  a  qTq.  of  full  rank,  and  thus  (5-1)  cannot  be 

J  J 


T""*;!!/  solved  for  the  elements  of  1^  for 

3 


positions  selected.  The 
Cf  this  chapter  is  to  prove  this  claim,  which  is  equivalent  to 


;  that  there  exists  an  infinity  of  solutions  for  the  elements  of  I.'/r 

J/  J 
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h  (5-1).  Note  that  (5-1)  is  based  solely  on  whole-body  mass  properties  for 
selected  positions  and  estimates  for  component  masses. 

The  method  of  proof  will  be  to  show  that  for  any  positions  I  and  II 
[only  j  different),  \l\  ^  has  the  same  linear  dependence  of  columns,  i.e. 
t.*e  sum  of  the  first  three  columns  is  zero.  Once  this  is  shown,  then  it  is 
trdly  seen  that  Q.  in  (5-1)  will  have  this  linear  dependence  since  each 

J 

the  submatrices  must. 

Consider  any  orthogonal  transformation  T,  which  for  this  work  is 
H  3.  By  the  definition  of  orthogonal,  =  T  ^ ,  and  consequently 
*•'  =  I.  Writing  this  out  in  terms  of  elements  and  performing  the  matrix 
implication: 


T(1 ,1 )  T(1  >2)  T(1 ,3) 
nT  =  T(2 ,1 )  T(2 ,2)  T(2 ,3) 

J ( 3 , 1 )  T(3,2)  T(3 ,3) 

i!:i'  lilently: 


T(1,D  T(2,l)  T(3,l ) 
1(1 ,2)  T( 2 ,2)  T(3 ,2) 

TO, 3)  T(2 ,3)  T(3 ,3) 


% 


|  p2(l»l)  j  [T(1,1)T(2,1)  j  [T(1,1)T(3,1)  | 

|  +T20,2)  +T(1 ,2)T(2 ,2)  +T(1 ,2)T(3 ,2) 

|  l+T2(l,3)j  ■  |+T(1 ,3)T(2 ,3)J  {+T(1 ,3)T(3,3)J 

i 

j  (T(2,1)T(1,1)  |  fT2(2,l)  I  fT(2,l)T(3,l)  |  [l  0  0 

j  !*T(2,2)T(1,2)  +T£(2 ,2)  +T(2,2)T{3,2)  =010 

j  (*T(2,3)T(1,3)J  l+T2(2,3)J  [+T(2,3)T(3,3)j  lo  0  1 

!  f’(3,l}T(l  ,1)  |  fT(3,l)T(2,l)  |  CT2(3,1)  1 

i  j*T(3,2)T(1.2)  +T(3,2)T(2,2)  +T2{3,2) 

f  ;‘T(3,3)T(1,3)J  |+T(3,3)T(2,3)J  +T2(3,3) 

for  any  positions  I  and  II: 


T‘(l,l>‘-+T‘(l,2)Sr(l,3)‘  .  T‘1(l,l)‘:+TJ10,2)<f+Tn(l,3)2  = 
TI(2,1)2+TI(2,2)2+TI(2,3)2  =  Tn(2,l)2+Tn(2,2)2+Tn(2,3)2  =  ' 
Tr(3,l)2+TI(3,2)2+TI(3,3)2  =  Tn(3,l)2+Tn(3,2)2+Tn(3,3)2  =  ' 

TI(2,1)TI(3,1)+TI(2,2)TI(3,2)+TI(2,3)TI(3,3) 

=  Tri(2,mn(3,l}+TII(2,2)TII(3,2)+Tn(2,3)TII(3,3)  -  t 

■TI(1,1)TI(3,1)+TI(1,2)TI(3,2)+TI(1,3)TI(3,3) 

'  TnO,l)Tn(3,l)+TnO,2)TII(3,2)+TII(l>3)Tn(3>3)  =  0 

TI(l,l)TI(2,l)+TI(l,2)TI(2,2)+TIn,3)TI(2>3) 

'  TI1(l,l)TII(2,l)+Tn0.2)T!I(2,2)+TII(l,3)TII(2,3)  =  o 

But  manipulating  the  above  equations,  and  then  substituting  the 
ents  of  V*  11  (referring  to  Equation  (4-8): 

V]"11(l,l)+Vj"11(l,2)+V]'11(l,3)  =  0 
Vj“II(2,l)+vj"II(2,2)+v]'II(2,3)  =  0 

Vj~  n(3’l  )+Vj_II(3,2)+Vj~rr  (3,3)  =  0 

vj‘II(4,l)+v]'II(4,2)+vj'II(4,3)  =  0 
yJI“II(5,l)+v]‘II(5,2)+v]“II(5f3)  =  0 
v j” 1 1 (6 » 1 )+V 1 1 ( 6 , 2 )+V 1 1 (6 , 3 )  =  0 


It  can  be  seen  from  these  equations  that  the  sum  of  the  first  three 


J  of  V.  is  zero.  Thus  vj  11  is  of  rank  less  than  6  for  any 


•  -s  I  and  II.  The  choices  of  positions  I  and  II  are  arbitrary 


MI  ..I-III  ..I-N 


VJ 


V. 

J 


will  all  have  this  same  linear  depen- 
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Tv1"11 

force  of  columns.  Clearly,  then,  Q.  =  ! 

J 

VI-N 

L  o 

'I  ’.ince  the  sum  of  its  first  three  rows  is  zero,  and  thus  it  cannot  have  six 

■;*.tarly  independent  columns.  Hence,  for  any  number  of  positions,  (5-1)  will 

V,  yield  a  unique  solution  for  the  elements  of  I.)  .  In  order  to  obtain 

» -iif.ique  value,  then,  it  is  necessary  to  introduce  an  additional  equation. 

fii  approximation  of  the  sum  of  the  component  moments  of  inertia  can  provide 

independent  equation,  which  together  with  a  Q.  matrix  formed  from 

J 

.mf.cious  selection  of  three  body  positions  determine  a  unique  value  for  the 

,v . 

‘Verts  of  I  (as  was  shown  in  Example  3).  J:CI 

J/c,  K  « 

j  » 


CHAPTER  6 


Suggested  Methods  for  Estimation 
of  Component  Masses 

The  determinations  of  the  component  centers  of  mass  and  inertia  tensors 
Upended  on  estimates  for  the  masses  of  the  components.  The  purpose  of  this 
chapter  is  to  provide  two  possible  methods  for  making  this  estimate.  These 
ire  regression  equations  and  water  immersion  (with  an  assumed  density).  Both 
methods  depend  on  statistical  data  for  cadavers.  Inclusion  of  outer  garments 
suits  is  also  considered. 

ly  cress  ion  Equations 

Estimates  for  component  masses  can  be  obtained  from  regression  equations 
fi- dieting  component  weights  from  total  body  weight. 

J.  T.  Barter  (Ref.  1)  used  the  frozen  cadaver  data  of  earlier  research- 
mi  to  devise  a  system  of  regression  equations  estimating  the  component 
*i*'/.ts.  These  equations  are  as  follows,  where  all  weights  are  in  Ibf. 
i-'.  11,  page  171): 
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'*^6> 

;.\Vv. 


Body  Segment 

Regression 

Equation 

Standard  Deviation 
of  the  Residuals 

Head,  neck  and  trunk 

= 

0.47 

X 

Total 

body  wt. 

+ 

5.4 

(±2.9) 

Total  upper  extremities 

= 

0.13 

X 

Total 

body  wt. 

- 

1 .4 

(±1.0) 

Both  upper  arms 

= 

0.08 

X 

Total 

body  wt. 

- 

1 .3 

(±0.5) 

Forearms  plus  hands3 

= 

0.06 

X 

Total 

body  wt. 

- 

0.6 

(±0.5) 

Both  forearms3 

= 

0.04 

X 

Total 

body  wt. 

- 

0.2 

(±0.5) 

Both  hands 

= 

0.01 

X 

Total 

body  wt. 

+ 

0.3 

(±0.2) 

Total  lower  extremities 

= 

0.31 

X 

Total 

body  wt. 

+ 

1.2 

(±2.2) 

Both  upper  legs 

= 

0.18 

X 

Total 

body  wt. 

+ 

1.5 

(±1.6) 

Both  lower  legs  plus  feet 

,= 

0.13 

X 

Total 

body  wt. 

- 

0.2 

(±0.9) 

Both  lower  legs 

= 

0.11 

X 

Total 

body  wt. 

- 

0.9 

(±0.7) 

Both  feet 

= 

0.02 

X 

Total 

body  wt. 

+ 

0.7 

(±0.3) 

li  =  11,  all  others  N  =  12, 


The  sum  of  the  predicted  weights  of  the  above  equations  does  not  always 
equal  the  total  body  weight.  To  compensate  for  this,  the  difference  is 
vilculated,  then  distributed  proportionally  among  the  segments. 

Another  set  of  regression  equations,  also  based  on  frozen  cadaver  data, 
mr.  developed  by  Clauser  and  McConville  (Ref.  4).  They  are  as  follows, 

»N'-re  all  weights  are  in  Ibf: 


Head  Weight 
Trunk  Weight 
Upper  Arms  Weight 
Forearms  Weight 
Hands  Weight 
Upper  Legs  Weight 
Lo wer  Legs  Weight 
Feet  Weight 


.028  Body  Weight  (W)  +  6.354 

.552  W  -  6.417 

.059  W  +  .862 

.026  W  +  .85 

.009  W  +  .53 

.239  W  -  4.844 

.067  W  +  2.846 

.016  W  +  1 .826 


I 

yi 
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As  for  Barter's  equations,  the  predicted  weights  are  corrected  so  that 
the  difference  between  total  body  weight  and  the  sum  of  predicted  weights  is 
distributed  proportionally  among  the  segments.  Note  that  Clauser  and 
McConville's  equations  provide  an  estimate  for  head  weight  not  given  by 
Barter's . 

It  is  claimed  that  this  second  set  of  regression  equations  results  in 
a  substantial  redistribution  of  weight.  Which  of  the  two  sets  of  regression 
equations  is  more  accurate  is  as  yet  an  unanswered  question. 

{ 

'i 

i; 

f  Water  Immersion 

|  By  having  the  subject  in  question  alternately  immerse  his  body  parts 

|  into  water  tanks,  the  volumes  of  components  can  be  measured.  Internal  compon- 

|  c-nts '  volumes  can  be  found  by  irmiersing  each  internal  component  with  the 

(>  corponent(s)  external  to  it,  then  subtracting  the  already  determined  volume(s) 

$ 

|  'f  the  external  component(s) .  The  volume  of  the  central  component  can  be 

\  found  by  immersing  the  entire  subject  to  find  his/her  whole-body  volume, 

If 

I  tv-n  subtracting  the  already  determined  volumes  for  the  other  body  components. 

I 

I  Using  statistical  data  for  specific  gravity  of  cadaver  parts  (Ref.  6, 

|  ' 

I  1955a),  approximate  values  for  component  masses  can  be  found  from 

§  <.t;*?riman tally  determined  volumes  by  assuming  that  the  subject's  components' 

ft 

|  -fcific  gravities  are  reasonably  close  to  the  cadaver  data.  However,  this 
f;  •.  not  account  for  individual  variations  in  body  structure.  Thus  if  one 
has  larger  than  average  bones  then  his  component  specific  gravity 
tv?  V;  expected  to  differ  from  the  average  data  for  specific  gravity. 

''  *'■'  method  of  water  immersion  may  be  qualified,  in  that  component 

ray  be  changed  by  water  pressures  encountered.  This  effect,  if 
;iiMe,  can  be  minimized  by  immersing  the  components  as  much  as  possible 
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near  the  surface  of  the  tank  (more  horizontal  than  vertical).  One  possible 
way  to  avoid  such  a  deformation  effect  would  be  to  use  rigid  casts  of  the 
subject's  body  parts  for  the  immersion  tests. 

This  water  immersion  technique,  though  more  difficult  and  costly  than 
using  regression  equations,  has  the  advantage  that  it  better  accounts  for 
individual  differences  in  mass  distribution. 

Inclusion  of  Outer  Garments  or  Sui ts  in  Component  Masses 

Each  component  mass  can  be  considered  to  include  the  mass  of  the  outer 
covering  surrounding  the  component  (see  Chapter  4).  The  contribution  to  the 
component  masses  by  the  outer  covering  parts  can  be  found  by  simply  cutting 
the  outer  covering  into  appropriate  sections  and  weighing.  These  values  are 
then  simply  added  to  the  corresponding  component  masses. 

If  destroying  a  garment  or  suit  by  cutting  it  into  sections  is  econom¬ 
ically  undesirable  or  impossible,  then  the  researcher  must  devise  his  own 
tothod  for  obtaining  reasonable  values  for  the  masses  of  the  outer  covering 
components . 

Tieber  and  Lindemuth  provide  a  set  of  regression  equations  for  estimat 
the  component  masses  for  a  typical  pressure  suit  (Ref.  17,  page  28). 


CHAPTER  7 


Suggested  Methods  for  Estimation  of 
The  Sums  of  Component  Moments  of  Inertia 

The  determination  of  the  component  inertia  tensors  depended  on  estimates 
for  the  sums  of  component  moments  of  inertia  (see  Chapters  4  and  5).  The 
purpose  of  this  chapter  is  to  provide  two  possible  methods  for  making  this 
estimate.  Methods  for  inclusion  of  outer  garments  or  suit  shall  also  be 
presented. 


Hanavan  Model 

The  first  method  is  to  obtain  values  for  s.  (j  =  l,...,n)  by  using 

J 

Hanavan's  simply  shaped  components  to  obtain  estimates  for  I yc  (j  =  l,...,n) 
(see  Chapter  3)  then  finding  s.  from  its  defining  equation: 

J 


;j  3  J,  V/c,  (k>  k) 


(j  1 . n) 


This  is  not  equivalent  to  assuming  simple  geometric  shapes  for  the 
body  components,  rather  it  is  assuming  that  the  sums  of  moments  given  by 
Hanavan's  model  will  not  differ  significantly  from  the  actual  values  for  the 


sums. 


Hanavan's  estimates  for  I.',  (j  =  l,...,n)  will,  as  was  indicated 

J  '  cj 

5r  Chapter  3,  depend  upon  measurements  of  the  subject's  dimensions  and 
estimates  for  component  masses.  The  researcher  can  use  either  of  the  methods 
resented  in  Chapter  6  to  obtain  component  mass  estimates. 
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Casts  of  the  Subject's  Components 

Molded  casts  of  the  particular  subject's  components  made  of  a  uniform 

density  material  provide  the  second  method  for  estimating  the  s  .  (j  =  l,...,n) 

3 

If  experimental  tests  are  performed  on  these  casts  identical  in  nature  to 
the  method  for  whole-body  tests  discussed  in  Chapter  4,  then  estimates  for 

I.1,  (j  =  l,...,n),  and  thus  for  s.  (j  =  n)  can  be  obtained. 

J  /Cj  3 

It  is  not  necessary  that  the  uniform  density  material  used  to  make  the 

casts  have  the  same  specific  gravity  as  the  data  for  body  components.  This 

is  true  because,  for  example: 


I 


xx 


J 


f  9  9 

[yc  +  z^)dm  = 


2  2 

(y  +  z  )  pdv  =  p 


J 


(y2  +  z2)dv 


m  v  v 

for  the  casts.  Thus  the  ratio  of  specific  gravities  can  be  used  as  a  scaling 
factor  to  obtain  estimates  for  I.',  (j  =  l,...,n)  for  the  subject. 

J  /cj 

It  is  worth  nothing  that  the  orientation  of  the  axes  about  which  the 
spin  tests  are  made  is  not  critical  as  long  as  the  three  axes  are  mutually 
perpendicular.  This  is  because  the  inertia  tensor  transforms  by  a  similarity 
transformati on.  It  can  be  easily  proved  that  for  an  orthogonal  similarity 
transformation,  the  sum  of  the  main  diagonal  elements  is  invariant  under 
ir»y  transformation. 


Inclusion  of  Outer  Garments  or  Suits  in  s. 

- . . 3 

There  are  two  methods  for  the  inclusion  of  outer  garments  or  suits  in 
devalues  for  s  .  (j  =  1 , . ,.  ,n) . 

J 

The  first  is  to  use  the  Tieber-Lindemuth  model  for  each  suit  part 
-'■-f-.-sponsing  to  a  body  component  (see  Chapter  3).  This  method  can  5c  aui 
find  the  inertia  tensor  for  each  suit  part,  which  is  then  simply  : 
ii  tt-o  inertia  tensor  for  the  cor  tv;  ponding  body  ccr  p-ancot  fc-t- 4  t  *?>*■•?* 
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by  using  Hanavan's  model  or  cast  experimental  data). 

The  second  method  involves  actually  cutting  up  the  garment  or  suit  into 

segments  corresponding  to  the  body  parts  and  performing  experimental  tests  to 

obtain  inertia  tensors.  This  can  be  done  separately  and  adding  in  or  actually 

surrounding  the  corresponding  cast  to  obtain  I.'/r  (j  =  1 . n)  for  the 

J  /Lj 

combined  body  component  and  suit  segment.  For  flexible  outer  coverings,  such 
test  data  may  be  of  questionable  validity.  For  a  pressure  suit,  it  may  be 
possible  to  inflate  each  component  to  sufficient  pressure  to  render  it 
sufficiently  rigid  to  solve  the  problem  of  flexibility.  This  would  involve 
somehow  sealing  off  the  ends  of  the  section  of  pressure  suit,  however. 

Conclusion 

This  chapter  has  presented  two  possible  methods  for  obtaining  reasonable 
estimates  for  the  sums  of  component  moments  of  inertia  (with  or  without 
outer  coverings  or  suits).  Both  methods  presented  depended  on  statistical 
data  for  cadavers  (Hanavan's  model  on  the  component  mass  estimates;  the  com¬ 
ponent  cast  method  on  the  specific  gravity  data). 

Hanavan's  model  is  clearly  the  easiest  and  least  expensive,  as  the  only 
experimental  data  needed  is  the  total  body  mass  and  the  various  anthropometric 
&sta.  But  the  component  cast  method,  while  more  costly  and  more  difficult 
to  carry  out,  offers  the  advantage  of  actually  taking  the  possible  variations 
from  statistical  average  of  component  shapes  into  account. 


CHAPTER  8 


Limitations 

There  are  four  types  of  limitations  on  the  method  for  predicting  mass 
properties  presented  in  Chapter  4.  These  are  errors  introduced  by  invalidity 
in  assumptions,  errors  in  whole-body  experimental  data,  accumulated  errors 
in  the  mass  property  predictions  for  the  central  component,  and  joint  con¬ 
straints  preventing  exact  application  of  the  method.  The  remainder  of  this 
chapter  discusses  these  and  presents  recommendations  for  further  study  where 
appropriate.  Finally,  suggestions  for  future  sensitivity  analyses  are  made. 

Invalidity  of  Assumptions 

There  is  some  question  as  to  the  validity  of  each  of  the  assumptions 
made  in  the  development  of  the  present  study's  method  for  determination  of 
component  mass  properties. 

As  was  indicated  in  Chapter  1,  the  rigid  body  assumption  is  of  marginal 
validity  for  the  torso  because  of  the  flexibility  of  the  backbone.  Hanavan's 
code!  partially  models  this  by  dividing  the  torso  into  two  halves.  Theoret¬ 
ically,  nothing  precludes  such  a  division  of  the  central  component  for  this 
study's  method.  It  remains  for  future  researchers  to  determine  the  validity 
Of  such  a  representation  of  the  central  component.  For  a  pressure-sui ted 
astronaut,  mobility  is  restricted  sufficiently  that  the  torso  may  prove  to 
it  adequately  modeled  by  only  one  component. 

The  hands  and  feet  are  also  flexible,  thus  also  compromising  the  rigid 
'*'■ '/  assumption.  But  these  body  parts  are  of  sufficiently  $rv>11  am  th*? 
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the  variation  in  contributions  to  the  mass  properties  due  to  flexibility 
may  be  negligible.  This  is  even  more  likely  for  a  pressure-suited  astronaut 
with  reduced  limb  mobility.  This  also  remains  for  future  researchers  to 
determine. 

The  assumption  that  there  exist  valid  methods  for  estimating  component 
masses  depends  on  there  being  little  significant  difference  between  component 
masses  of  a  live  subject  and  a  frozen  cadaver.  This  is  the  case  only  if 
component  properties  are  approximately  the  same  before  and  after  death  and 
if  individual  variations  are  not  too  great.  Some  indication  of  the  validity 
of  the  estimation  methods  would  be  given  by  comparing  results  after  applica¬ 
tion  of  both  methods.  Also,  because  the  total  number  in  the  sample  of  ca¬ 
daver  data  is,  to  date,  relatively  small,  there  is  a  need  for  more  such  data 
to  increase  the  validity  of  the  regression  equations  and  specific  gravity 
data.  Further,  as  was  indicated  in  Chapter  6,  the  effect  of  individual 
variations  in  body  structure  and  of  water  pressure  on  component  volumes  in 
water  immersion  experiments  must  be  considered  and  investigated. 

Next,  the  validity  of  methods  of  obtaining  sums  of  component  moments 
of  inertia  is  considered.  Both  methods  presented  {the  Hanavan  model  and 
cast  data)  require  estimated  values  for  component  masses  which  may  themselves 
is  in  error.  Also  both  methods  use  uniform  density  components,  which  actual 
iusan  components  are  not.  That  there  is  negligible  difference  is  another 
'*.tject  for  future  study. 

Finally,  the  assumption  that  the  components  of  the  outer  covering  are 
flexible  and  that  there  is  negligible  relative  motion  between  each  of  the 

components  and  the  corresponding  outer  covering  component  introduce 
r-]ikible  error  should  also  be  experimentally  investigated. 
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Errors  in  Whole-Body  Experimental  Data 

Errors  in  whole-body  experimental  data  arise  from  shifting  body  fluids, 
experimental  measurement  errors  (in  oscillation  frequencies  or  angular  rates 
and  torques),  movement  by  the  subject,  and  the  environment  of  the  experi¬ 
ments. 

The  error  due  to  shifting  body  fluids  can  be  minimi  zed, . as  was  indicated 
in  Chapter  4,  by  small  oscillation  amplitudes  and  long  pendulum  length  (for 
the  complex  pendulum)  and  by  small  angular  rates  (for  Fowler's  method). 

Future  work  would  involve  improvement  of  the  complex  pendulum  accuracy  and 
development  of  a  device  capable  of  physically  realizing  Fowler's  method. 

Inadvertent  movements  by  the  subject  can  be  minimized  by  careful  design 
of  the  apparatus  used.  Some  movement  is  unavoidable,  however,  specifically 
internal  movement  of  blood  and  expansion  and  contraction  of  the  lungs.  The 
influence  of  these  must  also  be  investigated  in  the  future. 

The  environment  of  the  whole-body  experiments  may  also  be  of  some  im¬ 
portance.  Probably  most  significant  in  this  respect  is  the  effect  of  gravity. 
If  the  experiments  are  carried  out  in  the  presence  of  appreciable  gravity 
(they  must  be  for  the  complex  pendulum),  then  the  gravitational  force  will 
tend  to  settle  the  body  fluids  toward  the  lower  portion  of  the  body.  Thus 
an  astronaut  under  weightless  conditions  may  have  an  appreciably  different 
mass  distribution  than  on  earth.  This  effect  may  be  negligible,  but  should 
be  investigated.  One  avenue  of  study  would  be  to  perform  the  experiments 
with  the  subject  both  upside  down  and  right  side  up  and  compare  the  resulting 
data  for  mass  properties. 

k cumulated  Errors  in  Central  Component 

As  was  indicated  in  Chapter  4,  the  central  component  mass  properties 
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can  be  determined  only  after  the  mass  properties  of  all  other  components 
have  been  found.  Thus  any  errors  made  in  the  determination  of  the  component 
mass  properties  will  accumulate  in  finding  the  mass  properties  for  the 
central  component.  There  is,  of  course,  the  possibility  that  errors  in  the 
component  mass  properties  will  cancel  out,  but  this  cannot  be  assumed. 

Joint  Constraints 

The  bending  of  certain  body  joints  is  constrained  so  as  to  preclude 
the  orientation  of  the  external  components  in  three  positions  not  in  one 
plane  (this  was  suspected  in  Chapter  4  as  necessary  to  uniquely  determine 
the  component  mass  properties).  The  elbows  and  knees  (hinge  joints)  are 
joints  of  this  type. 

For  the  joints  immediately  external  to  these  joints,  then,  a  researcher 
can  only  obtain  approximate  values  by  slightly  rotating  the  components  in¬ 
ternal  to  the  hinge  joint  for  one  of  the  selected  positions  so  as  to  allow 
the  external  component  to  get  out  of  the  plane  of  the  other  positions. 

Because  both  the  thighs  and  upper  arms  appear  fairly  symmetric,  this  seems  a 
reasonable  solution  to  the  problem  of  joint  constraints.  Experimental  data 
on  cadaver  thighs  and  upper  arms  can  serve  as  a  more  quantitative  measure  of 
this  symmmetry. 

Suggested  Sensitivity  Analyses 

Sensitivity  analyses  on  the  effect  of  the  possible  sources  of  error 
presented  in  this  chapter  should  be  carried  out  in  the  future  so  as  to 
determine  accuracy  requirements  for  reasonable  predicted  results.  These 
would  consist  of  alternately  determining  the  relationship  between  input 
errors  in  component  masses,  sums  of  component  moments,  and  whole-body  data 
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and  the  resulting  differences  in  predictions  output.  Of  secondary  importance, 
though,  is  the  differences  in  predicted  mass  propertires;  more  important  are 
the  differences  in  predicted  dynamics  of  the  system  of  interest. 


CHAPTER  9 


Conclusion 

In  this  study  a  method  has  been  presented  for  predicting  the  component 
mass  properties  for  a  human  about  component-fixed  axes,  with  or  without 
outer  coverings.  This  prediction  depended  on  a  component  rigid  body  assump¬ 
tion,  estimates  of  component  masses  and  sums  of  component  moments  of  inertia, 
and  data  for  whole-body  mass  properties  for  selected,  known  positions.  If 
an  outer  covering  was  to  be  included,  then  the  prediction  additionally 
depended  on  the  assumptions  that  negligible  error  was  introduced  by  consider¬ 
ing  the  covering  components  rigid  and  relative  movement  between  body  and 
covering  components  negligible. 

The  researcher  is  free  to  choose  the  number  and  definitions  of  the 
body  components.  Thus  if  it  turns  out  that  significant  error  does  not  re¬ 
sult,  the  hands  might  be  considered  part  of  the  forearm,  the  feet  part  of 
the  calves,  or  the  head  and  neck  part  of  the  trunk.  It  remains  for  future 
researchers  to  investigate  the  validity  of  different  component  definitions. 

Limitations  of  the  results  including  several  possible  sources  of 
error  were  discussed  in  the  previous  chapter.  Some  of  these  error  sources 
are  not  easily  isolated  nor  the  degree  of  error  easily  ascertained.  However, 
the  best  indication  of  the  validity  of  the  method  presented  in  this  study 
is  not  the  accuracy  of  the  predicted  mass  properties,  but  the  degree  of 
agreement  between  predicted  and  actual  dynamics  of  the  system  of  interest. 

In  summary  of  the  method  presented  in  this  study,  a  flow  chart  is 
provided  in  Figure  4  indicating  the  sequence  of  steps  taken  in  actually 
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carrying  out  a  prediction  of  the  mass  properties  for  a  component.  Dotted 
lines  indicated  optional  steps.  Once  component  mass  properties  have  been 
determined,  whole-body  mass  properties  can  be  predicted  for  any  position 
by  using  Equations  (2-7)  and  (2-13). 
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,  Determine  Outer  . 
Covering  Component  Mass1 


Figure  4.  Flow  Chart  of  the  Determination 
of  Component  Mass  Properties 


APPENDIX  A 


PROGRAM  WHOLE  (FORTRAN  IV) 


The  purpose  of  PROGRAM  WHOLE  is  to  calculate  whole-body  mass  properties 
from  component  mass  properties  about  component-fixed  axes.  An  explanation  of 
the  symbols  used  in  the  program  follows: 


Computer  Symbol 
A 


Study  Symbol  and/or  Explanation 
The  Point  A. 


FCTR 

Conversion  Factor  from  Degrees 

G 

The  Angle  y 

I  BA ( 3 , 3 ) 

The  Matrix  Ig^  . 

I C  A  ( 3  s  3 ) 

The  Matrix  1..^  . 

ICC (3,3) 

The  Matrix  I./r  . 

1/Li 

I CCP (3,3) 

The  Matrix  I . '  /r  . 

/Li 

The  Scalar  m^  . 

MB 

MI 

The  Scalar  m.  . 

P(3,3) 

The  Matrix  P.  . 

PI 

TT 

PRMAT 

A  Matrix -Printing  Subroutine. 

RBA(3) 

The  Vector  R^^  . 

RCA (3) 

The  Vector  R^  ^  . 

RCJ(3) 

The  Vector  R„  , ,  . 

_ci/Ji 

RC JP (3) 

The  Vector  R  '  . .  . 

_Ci  /Ji 

RJA(3) 

The  Vector  R,  . 

Ji 

The  Matrix  T.  . 

T  ( 3 , 3 ) 

TICCP(3,3) 

The  Matrix  T.  Times  1 . 1  /r  . 

1  1  /  v  2 
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U 

X 


The  Angle  y. 
The  Angle 


The  following  are  external  subroutines  (listing  is  not  provided  for 
brevity): 

External  Subroutine  Purpose 

MSAB  Matrix  Multiplication  (AB). 

MSABT  Matrix  Multiplication  (ABT). 

MATADD  Matrix  Addition. 

MSNULL  Sets  all  the  Elements  of  a  Matrix  to  zero. 


A  listing  of  PROGRAM  WHOLE  follows. 


60 


PROGRAM  WHDLF  (INPUT,  OUTPUT! 


«  THTS  PROGRAM  CALCULATES  the  FOLLOWING  * 

*  tt 

«  1)  COMPOMFNT  aNO  WHOl.F-ROnY  CENTfP  OF  MAGS  «• 

0  LOCATIONS  pEIATTVF  TO  aPPTTPAOY  POINT  a,  # 

•» 

?>  COMPoNFNT  AND  WHOLF-ROnY  IMFRTTA  MATRICES  •» 

*  ABOUT  THF  AXFS  SYSTEM  aT  A.  « 

*  « 

*  UNITS  -  « 

*  <> 

«  A  LOCATION  -  * 

«  •» 

*  ROOY  position  —  # 


DIMENSION  RC JP ( 0 ) ,  Pr J ( 3 )  ,  PCA(3)*  PJA<R),  RBA(3), 

P  TH»3)  t  P  ( 3 » 3 )  ,  TICPP(3,3) 

RFAL  mi,  mb,  T  ORp  (  3 ,  R )  .  T  C  C  ( 3  ♦  3 )  ,  TCA(3,3),  T  BA  (  3  *  3> 

PT  =  3* 141 5B?GSRSSR70 
FCTR  =  Pl/lflo. 

RF  A  0  POO,  MB,  w 

CALL  MSN1JI.L  (PRA,  3,  1.  3) 

CALL  MSNULL  (TPA,  3,  3*  3) 

PR  TNT  10? 

I  0?  FORMAT  ( 1 H 1  > 

no  m  *  u  n 

READ  ? 0 3 ,  MI 
PFAO  PC  1  ,  X,  G,  IJ 
PRTMT  100,  X,  G,  U 

X  =  X  <>  F  C  T  R  s  G  =  G*FCTR  $  U  =  U*FCTR 
100  FORMAT  (SX,  3F i 0,6) 

ST  NX  =  STN(X)  «  COSX  =  CoS ( X ) 

ST  MG  =  SIM(G)  F  COSG  =  CoS ( G ) 

STNU  =  SIN ( I ) )  t  COS"  =  CoS (U) 

T  (1  ,  1 1  =  COSXorOSG 

T  ( 1  ,  ?  )  =  COSX*STnG«STN!I  -  STMX^CoSU 
T  <  1  ,  3 >  =  COSX«-cTNG»roStl  +  STMXeSTNU 
T ( ? , 1 )  =  SiMx#rosG 

T  {•?,?)  =  STMX«STMG*StN|I  +  COFXttCoSU 

T ( R , 3 )  =  STWX«FImG#CRSU  -  COSXoSjNU 

Tn,l)  =  -SING 

T (  3  »  P  )  =  COSG*FTNU 

T  (3»  3)'  =  C0SG»OOSU 

READ  POP,  RCJP 

PFAO  ?0?,  RJA 


C  At.L  MS  AB  ( T ,  o  •  3  *  RCiP*  3*  1*  RCJ,  3*  3,  3) 
CALL  MATAnn  (POJ,  RJA,  PC A  *  3*  1*  3*  3*  3) 
no  7  K  =  1  ,  3 

RRA(K)  =  RRA(K)  *  PCA  (K)  #  (MT/MR) 


7  COMTINUF 


no  3  K  =  1  ,  3 

PF AO  nctP(K,L> *  L=1,r) 

3  CONTINUE 

call  MSAB  (T,  t.  3,  TCrP,  3,  3*  TICCP,  .3 »  3*  3) 
■  CALL  MSABT  (TIcCP*  3.  3,  j*  3*  3*  ICC,  3*  3,  3) 

P(l,l)  =  MJ  #  fDOA  (?) *PCA  (?)  *  RpA  (3)  ORrA  (  3)  ) 

P(P,?)  :  wi  #  (OCAm*PCAM>  *  PCA  (3)«PcA  (3)  ) 

P(3»3)  *  MI  «  (RCA(n*0CA(1)  ♦  PC*  (?)  «RcA  (?)  ) 

P<1  ,?)  =  MI  *  fPrAd  \*RCA<2)  ) 

P('l,3)  =  mt  *  (RCA  ( 1)  *PCA  (  3)  ) 

P<p,3)  S  MI  *  (RCA  (?) *PCA (3) ) 

P  <?,I  )  =  P (1 ,?) 

P  (  3  ,  1 )  =  p  M  ,  3 ) 

P(B,?)  =  p  (  ? ,  3 ) 

CALL  MATAnn  (IrC,  P,  ICA,  3,  3,  3*  3*  3) 

CALL  MATAnn  ( I  n  A  «  I C  A ,  IRA*  3*  3,  3*  3*  3) 

\  CONTINUE 


cai  l  prmat  < rr a »  3. 

CAI  I.  PRMAT  ( IRA  .  3, 


?oo  format 
?0I  FORMAT 

pop  format 

?03  FORMAT 
FNO 


(F)0. 6,  T 1 0  ) 
( 3  F  )  0 .  ? ) 
(3E?0.6) 

(FI  O.ft) 


i ) 
3) 


SURROIITINF  PRMAT  UTRlX*  K»l.) 

«  this  subroutine  prtmtr  any  k  rv  l  matrix.  # 

n t mfnston  atriy  (k,u 
PRTMT  101 

no  1 1  m  =  1  ,k 

PRTMT  100,  (ATpTX(m,m),  n  =  I,L) 

I  CONTINUE 

100  FORMAT  dx,  f>r?0.1?) 

I  0 1  FORMAT  (1/,  IX,  71  (IH«)  ) 

RFTtIRN  %  FMO 


APPENDIX  B 


PROGRAM  PART  (FORTRAN  IV) 


The  purpose  of  PROGRAM  PART  is  to  calculate  the  mass  properties  for 
any  body  component  (except  the  central  component)  from  whole-body  mass  property 
data  for  three  selected  positions,  an  estimate  of  the  component  mass,  and  an 
estimate  of  the  sum  of  the  component  moments  of  inertia.  An  explanation  of 
important  symbols  different  from  or  in  addition  to  those  symbols  of  PROGRAM 
WHOLE  (Appendix  A)  follows: 


Computer  Symbol 


Study  Symbol  and/or  explanation 


A(6 ,3) 

The  Matrix  A. 

J 

ATA(3,3) 

The  Matrix  A Ja. 

J  J 

ATAAT(3,6) 

The  Matrix 

1 

J  Jj 

■v 

BBI (6) 

A  Dummy  Vector  Needed  by  Subroutine 

MATINV. 

BBR(3) 

A  Dummy  Vector  Needed  by  Subroutine 

MATINV • 

BI (1 3) 

i 

The  Vector 

1 

b—t 

1— H 

1 

1 — 1 

1^ 

» _ 

T~  I -I  I  I 

bj 


Sj 


BR(6) 


The  Vector 


A I  —  1 1  /  \ 

a  (Rg) 


DETER 


The  Determinant  of  the  Matrix 
Inverted  by  MATINV. 
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The  Prefix  E 

The  Computer  Symbol  Prefixed  is  for 
a  Component  External  to  the  Component 
Whose  Mass  Properties  are  Sought. 

E(6)  (=E1+E2) 

A  Vector  Used  for  Notational  Brevity. 

EO (6,3) 

A  Matrix  Used  for  Notational  Brevity. 

El (6)  (= (EO) (ERCJP) 

A  Vector  Used  for  Notational  Brevity- 

E2(6) 

A  Vector  Used  for  Notational  Brevity. 

EI(13) 

A  Vector  Used  for  Notational  Brevity- 

MC 

m . 

J 

The  Suffix  N 

The  Symbol  Suffixed  if  for  Body 

Position  N. 

NE 

The  Number  of  Components  External 
to  the  Component  Whose  Mass  Properties 
Being  Sought. 

Q(13,6) 

* 

QTQ(6 ,6) 

*T  *  ,  *T  *. -1 

V  0R  <VV 

QTQQT 

<Q*Tq*rV 

J  J  J 

SC 

so 

V12(6,6) 

V  I"11 

J 

VI 3(6 ,6) 

V  i-m 

J 
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The  following  are  additional  external  subroutines  (listing  is  not 
provided  for  brevity): 

External  Subroutine 
MATINV 
MSATB 


Purpose 

Matrix  Inversion 

T 

Matrix  Multiplication  (A  B) 


VSUBC 


Vector  Subtraction 
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PROGRAM  PART  ( T  mput  ,  OiiTpyT) 


«•  O 

O  THTS  PROGRAM  r>r TfRm T k,e  s  THF  MASS  PROPERTIES  FOR  o 
O  ANY  COMPONENT  OF  TRF  ROD Y  EXCEPT  FOR  THf  CENTRAL  o 
o  COMPONENT  (TYPTOALI.Y  THE  TORSO).  0 
o 


DIMENSION  a  <ft,-5>,  ATA<3*3)*  ATAAT  <3«6)  ,  -RBI  (6)  *  RBR(3)» 

P  RTM3)*  RR(ft).  F(M.  F0(ft,3)*  El(6)*  E?(ft),  EH13)* 

3  ETCCP(3.3)  *  fd(3,3,3),  FPC  JP  ( 3 )  *  ERJA<3*3)«  ET(3*3*3)* 

4  FTTTM(3.3).  FTTTTn.3.3),  fTN(3,3),  IBA(3*3,3), 

5  ICCP(3*3)*  T  MDFR  (  3  ,  ■? )  .  InDFT(6»p)»  IPIvR(ft)*  IPT  VI  (ft)  * 

6  p  ( 3  » 3  ♦  3 )  .  0(13, ft),  oTo(ft,ft),  QToQT(ft*l3)*  P(3)* 

7  R  R  A  (  3 , 3  )  ,  RCJo(3),  dJA(3)*  SUMR(ft),  SUmI<13)*  T(3,3,3), 

8  V/ 1  ?  (  6  i  6  )  ,  V 1  3  ( ft  ,  6 )  ,  VcTRi(F) 

REAI  mr,  MC»  Mr,  Tra.  TCCP 

PI  =  3 • 1 4  I 5RpftG3G8R7R 
FCTR  =  PI/ipo. 

PR  T  NT  P  04 

RE  AO  pOO,  MR,  ME 

REAP)  Pn3 ,  mc 

no  i  n  =  i ,  3 

RF A n  ?P1,  X,  G,  U 

PR  TNT  ?01,  X,  G,  U 

X  =  XftFCTR  «  G  =  G^FCTR  S  U  =  U&FCTR 

1  CALL  FORMT  (t,  n,  X,  G,  U) 

DO  ?  K  =  1  ,  3 

DO  ?  L  =  1  *  3 

A  ( K  ,  L  )  =  TM,K,L)  -  T  ( ? ,  K « L  ) 

2  A  (  <  +  3  »  L )  =  T  (  1  ,  K  ,  L  )  -  T  (  3  ,  K  ,  |  ) 

CALL  mSATP  (a,  f,  t.  A,  ft,  3,  ATA*  ft*  ft,  3) 

CAI..L  M  a  T  T  n  V  (  ATA  ,3,RrR,  OfDETFR*  IpIVR*  INoER,  3,  ISCAR) 

CALL  MSART  (ATA,  3,  3*  A,  ft,  3,  ATAAT,  3*  ft*  3) 

00  3  N  =  1  *  3 
RE  A 0  P02,  R 
DO  3  K  =  1 ,  3 

3  R  R  A  (  N  ,  K  )  =  P  ( K ) 

00  ft  K  =  1  *  3 

PRfK)  =  (MR/MC)  *  ( RR  A (  1  ,  k  )  -  RR  A ( P  «  K )  ) 

E  rr  (  K  ♦  3 )  =  (MR/MC)  *  (  RR  A  (  1  ,  K  )  -  RB  A  (  3  *  K  )  ) 

IF  (NF.EQ.O)  GO  TO  13 
DO  IP  !=  I,  ME 
RFAO  p  03  *  ME 
DO  ft  K  =  1 , ft 
ft  SI|MR(K)  =  0,0 
OO  7  n  =  1,3 
REAn  ?01*  EX,  eg,  El] 
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E  X  =  FX^FCTP  <r  E  fi  =  FG^FCTr  $  EU  =  EU<>FCTR 
7  CALL  FORMT  (FT,  M»  EX*  EG ,  EU) 

CO  fl  K  =  1 ,  3 
DO  H  |.  =  1 ,  3 

E 0 ( K  ,  1. )  =  (MF/wr)  <*  (fT(i*k*L)  -  ET (?»«*L) ) 

P  E0(K+3,L)=  ( ME  /mo)  #  (FT(1»K,L)  -  ET(3,K»D) 

RF ah  ?0?,  ERCJp 

CALL  MFAB  (FO,  7,  EPCJP,  3*  1*  El*  6,  3*  3) 

0  0  Q  M  =  1  ,  3 
RFAO  ?D 2,  R 
00  Q  K  =  1  ,  3 
9  E  R  J  A  (  M  ,  K )  s  R  ( K ) 

DO  1 1  K  =  1 .  3 

E?(K)  =  (mF/mc)  #(EpJA(1»K)  -Fr J A ( ? » K ) ) 

11  E? ( K  +  3 )  =  (MF/mo)  *(ERJA(1,k)  -FRJA (3,K) ) 

CALL  VADOC  (El,  E?,  f,  6) 
call  VAOOO  (SIJmR,  f,  SilMR,  A) 
no  ?i  k  =  i,  n 

?1  SUM!  (K)  =  0.0 
DO  Id  K  s  I,  l 
RFAO  ?02,  R 

CO  1 fi  L  =  1 ,  3 
IF  EirCP(K*|.>  =  R(l  ) 

no  ?  o  M*  l ,  3 
DO  19  K  =  1,  3 
DO  19  l  =  1,  3 
19  ETM(K,I.  )  =  FT(m,kt*L) 

CALL  MSAB  (FTM,  7,  3,  FlCcP.  3,  3*  ETI,  3,  3,  3) 

CALL  MSABT  (ETT,  3,  7,  ETKJ »  7,  3,  ETITN,  7,  3*  3) 

DO  ?3  K  s  l,  3 
no  ? 3  L  =  1,  3 
?7  ETTTT(N,K,L1  =  FTTTNfK.L) 

?o  call  PAXTH  (EP,  M*  FT,  ME,  FRCJP,  ERJA) 

DO  P2  k  =  l  ,  3 
K  P  s  4~K 

CALL  POR(KP,  Ko»  KP) 

E  T  ( K )  =  F  T  T  T  T (1  ,  K  ,  K  )  -  FTlTT (  ?  ,  K  ,  K  )  +  EP(1*K,K)  -  EP(?*K,K) 

El (K  +  ft) =FTITT ( 1 ,k*K)  -  ETITT(?,K,K)  ♦  EP ( 1 , K , K )  »  EP(3,K.K) 

??  ET  (X  +  Q)  =FTTTT  (  i  ,KQ,KR) -FTtTT  (3»KQ,KR)  ♦  EP  ( ).  ,  KO  ,  KR  ) -F  P  (  3  ,  KQ  ,  KP ) 
ET  (1 3)  =  n.n 

CALL  VADOC  (SIJmT,  FT,  FUM  T  »  13) 

1?  CONTINUE 

CAIJ.  VFURC  <RR.  SUMP,  nP,  ft) 

CALL  VSURC  (RT,  FUmi,  nj,  13) 

17  'CONTINUE 

CALL  mfah  (ATAaT,  3,  6,RR,  a,  1,  RCJP,  7,  6,  3) 

CALL  PPMAT  (PC  IP,  7,  1) 

CAM.  FORmv  (VI  7,  T,  )  ,  ?) 

CALI  FORM V  (VlT,  T,  i,  3) 

CALL  FORMO  (0,  VI?,  Vl7) 
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Cam...  pr'mat  (q,  n.  f,) 

CALL  MSATR  (Q,  13,  ft,  n*  1  3,  6»  oTO,  13,  13»  6) 

CALL  PPMAT  (QTo,  ft,  ft) 

CALL  MATU'V  (0TO,ft,RRT,n,nETFP,lPlVIf  iNnEI  »ft,  TSCAT) 

CALL  PPMaT  (QTo,  ft,  O 

CALL  MSApT  (QTo,  6,  a,  Q,  13,  6»  QTQQT,  6,  13»  ft) 

B  F  a  0  ?n?,  rja 

CALL  PPMAT  (PJa,  3,  1) 

PPTMT  ?03,  Mf 

Co  1 A  m  =  1  ,  3 

Call  -AXTr-  (P,  V;,  T,  vr,  PCJR*  PjA) 

CO  14  K  :  1,  3 

RF  AO  ?0j>,  R 
CO  14  L  =  1,  3 
1  A  I  R  A  (  L  ,  K  ,  |_  )  =  O  (|  ) 

DO  IS  K  =  1,  3 

PTfK)  =  JRA(l,^,^)  -  I»A  (  ?  ♦  K  ,  K  )  -  P(1*K,K)  +  P  (  ?  ,  K  *  K  ) 

RT (K+ft)=IRA(),«»K)  -  !RA(3«K,K)  -  P(1*K,K)  ♦  P(3*K*K) 

KP  =  4-K 

CALL  POP  (KP,  ko,  KO) 

Rl (K+3)=JRA() ,kO.KP)-IRA(?«ko»KP)“P(1,KO»KP)+P(?»KO«KR) 
IS  RI (K+P) =T»A ( 1 ,kO,KR)-IpA (3,K0,KR) -P ( 1 tKO»KR) +P (3»KO,kP) 
RFAO  ?A3,  SC 
RI ( 1 3)  =  SC 
call  ppmat  (RT ,  13,  1 ) 

call  WRAP  (OTQoT,  ft,  1 3 ♦  RI,  13»  1,  VCTPI,  6,  13«  6) 

DO  1 7  K  =  1  ,  3 
17  ICCP(K.K)  =  VCTPT(K) 

ICCP  (  1  1?-)  -  VC TPT  (4) 

ICOP ( 1 1 3)  =  VCTPT (si 
I  COP ( ?  »  3 )  =  VCTPT(ft) 

CALL  PPMat  (TCCP,  R,  31 


300 

FORMAT 

(FlO.fi, 

301 

FORMAT 

(3F10.P) 

30? 

FORMAT 

(3F?0,fi) 

303 

FORMAT 

(El O.ft) 

304 

FORMAT 

( 1 R 1  ) 

EMO 
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SIJRPOIJTIK'F  FORmT  ( T ,  m,  Xt  G*  II) 

#  THTS  SUBROUTINE  FORMS  TwF  MftTRTX  T.  ■» 


OTmFK'STOn  T (3,3*3) 

STMX  =  SIM(X)  s  COSX  =  CoS(X) 

STMfi  =  SIM(G)  t  COSO  =  CoS (0) 

s  i  mu  =  stm(U)  «  cnsu  =  cosdj) 

T (M, 1 , 1 ) =C0SX«r OSG 

T  (W,  1  ,?)  =  COSX*STMG*STNI)  -  STNX*C0SU 
T  (m«  1  «  3)  =COSXtfsTMG*COSU  ♦  STNiX^StNU 
T (M,?,l ) =siNx*rosG 
T  (m<?i ?)  =simx*stmg«stnm  ♦  cnsx#CoSU 
T  ( M * ? « 3 ) sST NX^s T mG#C0SU  -  COSX#SlNU 
T (M, 3,] ) =-SlMG 
T sCnS6#STMU 
T (Mi  3 . ?) =rosG#rosu 
RETURN  *  FmO 


SllRROUT  IMF  F0RM0  (0,  V)?,  V 1  3 ) 

*  This  S'JRROUTTMF  FORMS  TuE  maTRTX  o.  * 

4  (Httt9  999#<Ht«99«99*9«9«#99#9<HHHHHMHH}99ft 

DIMENSION  0(n,  0),  V 1  2  ( 6  *  6 )  *  V 1 3  (  6  *  6 ) 

Dn  i  k  s  I*  fi 


00  1  L  =  U  6 
Q  ( K  ♦  U  =  V 1  7  ( K  ,  I  ) 

G  |K+6,L)=V13  (K,l  ) 
QlHfl)  =  1 .  n  * 

0(13*4)  =  0,0  s 

Rfturn  S  f 


O ( 1 3,?)  =  l.o 
0  ( 1  3  *  G )  3  0.0 


0(13*3)  =  1.0 
0(13*6)  =  0.0 


StJRROUTlNF  FORmV  (V,  T.  M,  M) 

*  THIS  S'  IF?  ROUT  T  NF  FORMS  T  *<E  MATRIX  V.  * 

DIMENSION  V (6*G) ♦  T(0*0*3) 

00  1  K  =  1  ,  3 
DO  1  L  =  1  *  3 

1  V(x-.L)  =  T  (M,k,I  )  *T  (m*K*L)  -  T(N,K*L)*T(N*K,Ll 
DO  R  L  =  4*  6 
IP  =  7-L 

C At  L  POR  (LP,  |0.  |  R s 
DO  2  K  =  1  *  3 

?  V  (K*l)  =2.0*  (T(m,k*|.0)  *T  (M,K,|_R)  -T  <N,K,Ln)  *.T  (N,K,LP)  ) 

D  0  3  K  =  4  *  6 


KP  =  K-6 
CALL  PQR 
DO  7  L  = 

3  V (KiL)  = 
DO  4  K  = 
KP  s  K-3 
CA!  t  ROR 
DO  4  L  = 
LP  =  7-L 
CA!  L  DOR 

4  V  ( K  ,  L 1  = 
? 

rfturn 


(KP,  KO,  KR  > 

1,3 

T  (M,KO,l  )  *T  (M,KR,L)  -  T(N,KQ*U#T(N*KR,L) 

4 ,  f)  ■ 

(KP,  kD,  KRj 
A,  6 

( Lp ,  1.0,  l_R1 

T  (M,ko*|_R»«T  (M.KR.LO)  ♦  T  (M,KO,lO)  *T  (M*KR*LR) 
T  (N,Kn,|_Rl  #T  (m,kr.|_0)  -  T  ( N  *  KO  ,  LO )  ( N , K  R  •  LR ) 

*  EMIT 


SUBROUTINE  PAXTH  (p,  N.  T,  mc»  PCJP,  RJa) 

*  THTS  S'  IP  ROU  T  T  NF  ApD(_lr<;  THE  PARALLEL  AXIS  THEOREM «  * 

<HHt  *«#*'»•»*#***#  tt  *  *  tt  *  W  #*******#•»**#*  **##**•«•***##***#**## 

CTMFNSION  P(3,3,7),  RN(3,7)«  RCA{3>,  RCj(3>»  RCJP(3), 
?  RJA (3)  ,  T (3,3,3)  ,  T  M ( 7  »  3 ) 

R F  A 1.  wr 

PRTNT  1(11,  mc 

C  AH-  RWMAT  (PC  P,  7,  1) 

CAIJ.  PRMAT  ( RJ  a  ,  3,  1  ) 

DOT  K  =  1  ,  3 
DO  1  L  =  1  ,  3 
1  TM(K,L>  =  T  (  N  ,  x” ,  (  ) 

CAI  I.  PRMAT  (TN,  7,  71 

call  MSA9  (  TM  •  7 ,  3,  RcJP»  3*  1»  RCJ,  7'  3»  3) 

CALL  PRMAT  (RCJ.  3,  i ) 

call  vadoc  (rcj.  rja,  rca,  3) 

CALL  PRMAT  (RCA.  3.  1) 

no  ?  TP  =  1  .  3 
CALL  POP  (IP.  TO,  I R  1 
PRTMT  10n,  IP,  TO,  TR 

100  FODMAT  dox,  3  T  A 1 

P(M.TP.IP)  =  Mr*  (RC.A  (  TO)  *RCA  (  TO)  ♦  RCA.<  I'R)  *RCA(IR>  ) 
P(M.TO.TR)  s  Mr  *  (  Rr A ( TO) *RCA ( I R)  ) 

P (M. TR, IQ)  s  P (M, TO, tR) 

PRINT  01,  P ( N . T  P , T  P 1  ♦  P ( N , T  0 1 1 R ) 

101  FORMAT  (inx,  ?F?n.i?i 
?  CONTINUE 

DO  3  K  =  1 ,  3 
DO  7  L  =  1  ,  3 
7  PM ( K , L )  =  P ( N , K ,  |  ) 

CAI  I.  PRMAT  (PN,  7,  31 
RFTl'RM  %  END 
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StJ«BOUTTMF  PQR  (TP.  TO,  IR) 

*  This  RUpROUTTNF  nc-TFRMTMEF  THE  INTEGERS  Q  and  p  from  p.o 

IF  (TP.EQ.l)  Gn  TO  1 

IF  (TP. EG.?)  GO  TO  ? 

IF  (JP.EQ.R)  GO  TO  7 

I  1 0  =  ?  «  TR  =  -5  S  GO  TO  4 

?  TO  =  I  %  T  R  =  0  S  00  TO  4 

T  10=  1  T  IR  a  7 

4  CONTINUE 

RrTURN  T  FMO 


S'IIPPOUTImF  PRMaT  (ATOIV,  KtU 

(MHHHUHMl  O  #  -tv  «  ■»  «•  <HH>  <j.  #  «•  #  •&  #  <*<><)■  •&<!■«  *  #  ttfl.  •»  Vt 

«■  this  currouttne  potmTr  any  k  ry  l  matrix,  # 

OIMFNJSTON  ATRTX  (K,l_) 

POINT  101 

DM  M  s  1»K 

1  PR TNT  10n,  { A  T  R  T  X ( w , M ) ,  N  =  i»L) 

100  FORMAT  (IX,  *Sr?n.l?) 

1  01  FORMAT  (]  /,  IX,  71  nu*)  ) 

RFTIJRM  <E  END 
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